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ABSTRACT 


Terahertz  (THz)  radiation  is  a  portion  of  the  eleetromagnetie  speetrum  with  many 
potential  uses  but  without  a  wide  availability  of  powerful  sourees.  The  free  eleetron  laser 
(PEL)  has  been  demonstrated  to  produee  intense,  eoherent  THz  radiation.  This 
dissertation  explores  a  novel  THz  FEE  oscillator  design  that  is  tunable  within  a 
wavelength  range  of  ~60  micrometers  to  '-ISO  micrometers  and  could  produce  up  to 
kW  average  output  power.  It  utilizes  superconducting  spoke  cavities  for  the  accelerating 
structures  to  generate  an  8  MeV  electron  beam  with  ~8.75mA  average  current.  The 
variable  gap  undulator  has  only  10  periods;  this  design  choice  enhances  the  extraction 
and  therefore  output  power,  but  also  reduces  single  pass  gain.  To  overcome  the  reduced 
gain,  the  optical  cavity  is  designed  to  minimize  round  trip  losses.  This  design  is  relatively 
compact  and  could  be  installed  on  a  ship  to  test  technologies  that  are  relevant  for  future 
scaling  to  weapons  class  output  powers. 
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I.  INTRODUCTION 


A.  MOTIVATION 

Free  eleetron  lasers  (FELs)  are  eontinuously  tunable  sourees  of  eoherent 
eleetromagnetic  (EM)  radiation.  They  ean  be  designed  for  a  wide  range  of  industrial, 
aeademic,  medieal,  and  military/homeland  seeurity  applieations.  This  dissertation  will 
touch  on  the  applications  of  terahertz  (THz)  radiation  and  be  centered  on  the  design  of  an 
EEE  in  the  THz  regime.  THz  radiation  covers  a  frequency  range  from  -'300  Gigahertz 
(GHz)  to  10  THz  (wavelengths  of  -1  mm  to  300  pm).  This  particular  frequency  band 
covers  what  is  known  as  the  “THz  gap”  where  there  are  few  coherent  sources  and 
detectors  as  compared  to  other  widely  used  frequency  bands  both  above  (microwave  / 
radio  waves)  and  below  (infrared)  [1].  This  gap  exists  in  part  due  to  the  scarcity  of  viable 
conventional  gain  media  with  energy  transitions  that  correspond  to  THz  frequencies.  This 
dissertation  will  include  discussions  of  the  uses  for  THz  radiation,  comparisons  between 
the  different  sources,  simulations  of  the  gain  process  for  a  THz  EEE,  and  a  plausible 
physical  design  for  a  THz  EEE. 

B.  USES  FOR  THZ  RADIATION 

THz  radiation  has  a  wide  range  of  applications  for  several  professional 
communities. 

I.  Academic  /  Scientific  Applications 

Scientific  applications  include  THz  spectroscopy  and  pump-probe  experiments. 

a.  THz  Spectroscopy 

Molecular  spectroscopy  can  utilize  THz  radiation  to  study  various  materials  due 
to  the  strong  THz  absorption  lines  corresponding  to  molecular  vibrational  and  rotational 
modes  [2].  The  most  common  method  for  these  measurements  uses  relatively  weak 
incoherent  black  body  sources  and  Eourier  transform  infrared  spectrometers  (ETIR)  as 
part  of  the  detectors.  New  methods  may  be  developed  using  THz  EEE  to  produce 
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powerful,  coherent  sources  to  study  and  categorize  specific  spectral  line  emissions  from 
single  vibrational  or  rotational  modes. 

b.  Pump-probe  Experiments 

Pump-probe  experiments  can  be  utilized  to  measure  ultra-fast  phenomena  in 
materials.  The  general  principal  is  that  a  “pump”  light  pulse  (usually  a  high  peak  power 
laser)  is  incident  upon  a  sample  and  causes  an  excitation.  A  certain  time  later  a  “probe” 
pulse  (generally  lower  in  power)  is  incident  upon  the  sample,  and  the  reflection  or  the 
transmission  signal  is  then  be  measured.  A  current  area  of  research  at  several  national 
labs  is  to  use  a  THz  pump  source  along  with  an  x-ray  probe  source  [3].  The  timing 
between  the  two  pulses  is  altered  so  that  the  transient  behavior  of  the  relaxation  can  be 
measured  [4].  The  National  High  Magnetic  Field  Laboratory  is  interested  in  a  THz  FEL 
to  use  in  conjunction  with  high  peak  magnetic  fields  for  research  into  new  material 
properties.  Broadband  thermal  sources  commonly  employed  in  THz  spectroscopy  are  not 
suitable  for  pump-probe  experiments;  however,  the  THz  pulses  produced  by  an  FEL 
could  be  utilized  for  these  kinds  of  experiments. 

2.  Homeland  Defense  Applications 

A  useful  property  of  THz  radiation  that  lends  itself  to  many  uses  in  security  and 
homeland  defense  applications  is  that  it  passes  through  many  materials  that  are  opaque  to 
radiation  in  the  optical  regime.  THz  can  readily  penetrate  most  non-metallic  and  non¬ 
polar  materials  [5].  These  include  cloth  (clothing,  bags),  building  materials  (paper, 
cardboard,  wood),  and  a  wide  range  of  other  barriers,  which  may  be  used  to  conceal  a 
potential  threat.  Photons  in  the  THz  regime  are  of  low  enough  energy  to  be  non-ionizing, 
so  that  the  EM  radiation  itself  is  a  minimal  health  risk  to  the  person  being  scanned  or  the 
operator  of  the  device  [6].  This  type  of  technology  is  already  utilized  in  current  airport 
security  scanners. 

The  THz  radiation  tends  to  be  highly  absorbed  by  water  and  O2  causing 
propagation  in  the  atmosphere  to  be  fairly  limited  (see  Eigure  1).  However,  a  more 
powerful  THz  source  that  could  be  tuned  to  a  frequency  of  lower  atmospheric  absorption 
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could  open  up  the  door  to  longer-range  applications  of  THz  radiation.  An  example  of 
atmospherie  attenuation  of  THz  radiation  ean  be  seen  in  Figure  1. 


Extinction  spectrum  of  THz  radiation 


Figure  1 .  Atmospheric  absorption  bands  created  using  LEEDR  atmospherie 

absorption  code. 

For  example,  a  ~1  kW  THz  souree  can  transit  several  hundred  meters  and  still 
maintain  a  high  enough  signal  to  be  measured  if  tuned  to  a  low  absorption  band.  THz 
sources  that  are  eurrently  used  tend  to  have  many  orders  of  magnitude  less  output  power 
and  are  generally  not  tunable. 

Another  seeurity  application  of  THz  radiation  is  to  actually  utilize  its  poor 
transmission  properties  in  the  atmosphere  to  make  a  secure  communieation  system.  The 
high  attenuation  of  the  THz  radiation  and  the  stronger  direetionality  (over  submillimeter 
or  radio  waves)  eould  allow  for  a  more  secure  connection  between  linked  deviees  [7]. 

Furthermore,  THz  speetroscopy  could  be  employed  in  the  search  for  explosive  or 
illegal  substanees.  Many  explosives  (TNT,  C-4,  HMX,  RDX,  ete.)  and  illegal  drugs  (e.g., 
methamphetamines)  have  a  eharaeteristic  absorption  speetra  in  the  THz  region  [8].  This 
oharaeteristie  eoupled  with  the  ability  to  penetrate  many  dielectrie  materials  eould  aid  in 
the  seareh  for  eoneealed  illegal  substanees. 
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C.  TERAHERTZ  SOURCES 

Several  different  methods  are  used  to  ereate  THz  radiation  of  various  power  levels 
and  frequeneies;  a  few  of  them  will  be  diseussed  in  this  seetion.  One  main  souree  of  THz 
radiation  is  the  infrared-pumped  gas  laser,  whieh  produces  between  1  mW  -  20  mW  of 
continuous  wave  (CW)  power.  This  technology  has  been  proven  in  the  500  GHz  to  3  THz 
regime  but  the  power  spectrum  depends  on  the  gas  being  used  as  the  gain  medium  [9]. 

Quantum  cascade  lasers  are  a  relatively  new  source  that  is  starting  to  become 
mature.  These  semiconductor  devices  can  supply  ~1  mW  to  8  mW  average  power  and 
can  be  pulsed  at  10  kHz.  They  have  the  ability  to  be  tuned  from  -'20  GHz  to  -200  GHz 
[10]. 

Optical  down-conversion  is  a  method  where  a  higher  frequency  signal  (optical)  is 
converted  to  the  THz  regime  through  certain  interactions.  One  particular  method  uses  two 
high  intensity  optical  pulses  focused  onto  a  pair  of  electrodes  with  a  gap  <  -  1  pm  printed 
on  a  semiconductor.  This  produces  a  photocurrent  that  is  modulated  at  the  frequency 
difference  between  the  two  incident  laser  pulses.  This  is  tunable  to  produce  a  THz  source 
with  average  power  levels  -0.1  pW  to  1  pW  [11].  Other  systems  use  a  non-linear  crystal 
to  mix  frequencies  of  optical  pulses  to  produce  THz  radiation,  but  often  these  have 
average  power  levels  -nW  to  pW  [2]. 

Electron  tube  sources  are  available  that  generate  THz  radiation  from  the  emission 
of  bunched  electrons  that  spiral  in  a  strong  magnetic  field.  These  can  take  the  shape  of 
backward-wave  electron  tubes  that  have  a  wide  frequency  tuning  range,  and  are  fairly 
compact  (table  top),  but  supply  only  moderate  (mW  level)  power  [12].  There  are  larger 
microwave  and  THz  sources  in  the  form  of  gyrotrons  that  can  be  designed  to  provide 
higher  power  (-100  kW  to  1  MW)  in  either  continuous  wave  or  pulsed  mode  [13],  though 
the  pulses  are  of  -1  ms  length  and  are  not  as  desirable  for  pump-probe  experiments. 
These  are  larger  room  size  installations  and  often  aren’t  scalable  throughout  the  THz 
spectrum. 
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Along  the  lines  of  these  sourees  are  FELs,  whieh  also  utilize  beams  of  relativistic 
free  electrons  to  generate  the  radiation,  and  can  be  scaled  up  to  relatively  high  power  (kW 
level)  with  shorter  pulse  lengths.  A  few  of  the  current  THz  FELs  are  listed  in  Table  1. 


FELs 

Average  Power 

Peak  Power 

Wavelength  Range 

Novosibirsk  (Russia) 

-500  W 

-0.5  MW 

120pm  -  230  pm 

Kaeri  (Korea) 

6  mW 

-I  kW 

100pm  -  1200  pm 

FELIX  (Netherlands) 

-I  W 

-100  MW 

2pm  -  250  pm 

ISIR  (Japan) 

-O.I  W 

-0.5  MW 

32pm  -  144  pm 

FELBE  (Germany) 

-40  W 

100  MW 

18pm  -  250  pm 

CLIO  (France) 

-I  W 

-100  MW 

3pm  -  150  pm 

Table  1.  Current  THz  FELs,  from  [14]-[20]. 


These  FELs  have  high  peak  power  and  average  power  relative  to  other  THz  sources. 

1.  Advantage  of  THz  FELs 

THz  FELs  have  many  advantages  compared  to  other  THz  sources.  They  can 
achieve  relatively  high  average  power  (-'1  kW-2  kW)  compared  to  other  THz  sources, 
have  narrow  linewidths,  are  linearly  polarized,  have  short  pulse  durations  of  several 
picoseconds,  produce  radiation  in  a  fundamental  transverse  mode,  and  are  tunable  over  a 
wide  range  of  wavelengths.  The  short  pulse  length  of  the  THz  pulse  along  with  the 
fundamental  transverse  mode  shape  allows  for  high  peak  brightness  when  the  beam  is 
focused.  The  short  pulses  and  high  repetition  rates  make  FELs  useful  for  pulse-probe 
experiments.  For  molecular  spectroscopy,  the  ability  to  tune  the  frequency  of  the  FEE 
would  allow  specific  excitation  of  specific  modes  to  be  studied  [21]. 

D.  TERAHERTZ  EEL  DESCRIPTION 

The  THz  FEE  described  in  this  dissertation  is  being  developed  by  the  U.S.  Navy 
(USN)  and  Niowave,  Inc.  The  designed  spectral  region  is  ~66  |um-114  qm.  The  FEE 
generates  light  by  transferring  energy  from  relativistic  electron  bunches  to  an  optical 
pulse.  It  does  this  by  transversely  “wiggling”  the  electrons  inside  a  periodic  magnetic 
field,  causing  the  electrons  to  radiate.  The  electron  source,  or  injector,  will  supply 
electrons  up  to  an  energy  of  -'100  keV.  Next,  the  electrons  will  be  injected  into  a  linear 
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accelerator  (linae)  strueture  that  will  bring  the  eleetron  beam  to  the  final  design  energy  of 
8  MeV.  This  relativistie  eleetron  beam  will  then  be  direeted  through  an  undulator 
(periodic  rows  of  alternating  permanent  magnets)  that  generates  the  periodie  transverse 
magnetic  field  that  eauses  the  eleetrons  to  radiate.  The  wavelength  emitted  depends  on 
the  electron  beam  energy,  the  strength  of  the  magnetie  field,  and  the  period  of  the 
magnetic  field.  This  PEL  will  be  '-lO  m-15  m  in  length  with  a  cross-sectional  area  of  '-I 
m  ;  it  will  operate  as  an  oseillator,  whieh  has  an  optieal  eavity  surrounding  the  undulator 
and  beam  line  allowing  the  fields  to  build  up  inside  the  eavity  over  many  passes. 

E.  NIOWAVE,  INC. 

Niowave,  Ine.,  is  a  eompany  in  Lansing,  MI,  that  speeializes  in  ereating  super- 
eondueting  eleetron  aeeelerating  struetures  and  supporting  teehnologies.  These  linacs 
aeeelerate  eleetrons  up  to  high  energies  (near  the  speed  of  light)  so  the  eleetron’ s  kinetie 
energy  ean  be  utilized  in  a  wide  a  range  of  applieations.  These  inelude  the  ereation  of 
medieal  radioisotopes;  the  production  of  light  in  free  eleetron  lasers,  x-ray  sourees,  and 
THz  sourees,  and  neutron  produetion.  Niowave  also  builds  radio  frequency  (RF)  guns, 
superconducting  cavities  and  eryomodules,  high-power  RF  couplers,  and  eommereial  4K 
refrigerators.  The  Navy  is  partieularly  interested  in  the  use  of  supereonducting  linacs  to 
drive  high  average  power  FELs  for  shipboard  defense  [22].  The  researeh  done  in  this 
dissertation  is  for  a  proposed  THz  FEE  built  by  Niowave  for  the  USN. 

F.  EEL  ON  A  TEST  SHIP 

This  relatively  small  THz  FEE  is  being  designed  for  eonsideration  to  be  put  on  a 
USN  test  ship,  the  USS  Paul  F.  Foster.  The  USN  is  interested  in  a  high  average  power 
FEE  for  taetieal  seenarios  sueh  as  ship  defense.  The  desired  FEE  would  be  '-I  MW 
output  power  at  a  wavelength  with  little  atmospherie  absorption  in  a  maritime 
environment.  The  FEE  being  designed  by  Niowave  and  NFS  has  a  mueh  lower  eleetron 
energy  leading  to  a  longer  wavelength  and  an  average  power  several  orders  of  magnitude 
below  the  final  design  goal  for  ship  based  weapon.  This  projeet  could,  however, 
demonstrate,  on  a  ship,  many  of  the  neeessary  technologies  for  a  weapon  class  FEE. 
Things  such  as  the  liquid  helium  eryoplants,  the  eryomodules  that  surround  the 
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accelerating  structures,  the  optical  cavities,  and  eleetron  beamlines  need  to  be  vetted  to 
show  that  the  movement  of  the  ship  will  not  prevent  the  devices  from  working  under 
various  sea-state  conditions.  An  PEL  is  sensitive  to  the  alignment  of  its  components, 
whieh  includes  the  optical  cavity,  the  eleetron  beam  and  undulator  axes,  as  well  as  the 
beam  line  eomponents  that  guide  the  electron  beam  through  the  FEE.  The  THz  FEE  will 
use  a  relatively  new  shape  for  the  electron  beam  aeeelerating  strueture,  the  spoke  eavity 
design,  that  deereases  the  overall  size  while  making  it  more  rigid,  whieh  should  be 
helpful  in  a  shipboard  environment.  The  low  eleetron  beam  energy  of  the  THz  FEE 
deereases  the  amount  of  shielding  that  is  required,  making  the  whole  system  smaller  and 
lighter  in  weight. 
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II.  BASIC  FEL  THEORY 


A.  FREE  ELECTRON  LASER  COMPONENTS 

An  FEL  oscillator  is  composed  of  four  major  components:  an  injector,  which 
generates  the  electron  bunches  at  a  low  beam  energy;  a  linear  accelerator,  which 
accelerates  the  electron  bunches  up  to  the  relativistic  energies  needed  for  the  FEL 
interaction;  an  undulator,  which  causes  the  electrons  to  emit  light;  and  an  optical  cavity, 
which  stores  the  optical  pulse  to  build  up  the  power.  A  schematic  of  an  EEL  oscillator 
can  be  seen  in  Eigure  2. 


Eigure  2.  Schematic  of  an  FEL  oscillator,  from  [23]. 


1.  Injector 

The  injector  is  made  up  of  two  key  components:  the  cathode  to  generate  the  free 
electrons  and  a  method  to  immediately  accelerate  the  low  energy  electrons  away  from  the 
cathode  surface. 

a.  Cathode 

An  EEL  requires  an  electron  source  (cathode)  that  will  produce  bunches  of  free 
electrons  with  proper  beam  characteristics  needed  for  the  EEL  interaction.  The  main 
types  of  cathodes  are  photoelectric,  field-emission,  and  thermionic.  Which  type  is  chosen 
depends  on  the  electron  beam  parameters  that  are  needed  for  the  EEL  interaction.  Eor  the 
proposed  Niowave  THz  FEL,  the  thermionic  cathode  will  be  implemented. 

A  thermionic  cathode  uses  the  thermal  energy  imparted  to  an  electron  in  a 
conduction  band  to  give  it  enough  energy  to  overcome  the  surface  work  function  of  the 
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material  and  become  unbound  and  enter  the  vacuum.  Thermionic  cathodes  have  several 
benefits  including  long  lifetimes  and  the  ability  to  supply  high  current.  A  couple  of 
disadvantages  of  thermionic  cathodes  are  that  they  generally  have  lower  beam  quality  and 
cannot  be  modulated  as  easily  as  other  types  [24]. 

b.  Accelerating  Cavity 

Space  charge,  the  inter-electron  Coulomb  repulsion,  can  cause  the  newly  created 
beam  to  rapidly  spread  out  in  space  as  well  as  energy.  If  left  unchecked,  this  can  destroy 
the  beam  quality.  However,  as  the  electron  beam  approaches  the  speed  of  light 
space  charge  effects  are  significantly  reduced  due  to  relativistic  effects.  Thus,  the  newly 
generated  electron  beam  should  be  quickly  accelerated  to  relativistic  energies  to  maintain 
good  beam  quality.  This  acceleration  is  done  using  either  a  radio  frequency  (RF)  or  direct 
current  (DC)  electric  field  at  the  cathode.  This  combination  of  cathode  and  accelerating 
cavity  makes  up  the  injector. 

A  DC  injector  uses  a  static  potential  between  the  cathode  and  an  anode  to 
accelerate  the  electron  beam.  A  major  problem  with  DC  injectors  is  that  the  constant, 
high  potential  gradients  can  cause  electrical  breakdowns  to  occur.  To  mitigate  this 
problem,  an  RF  cavity  is  often  utilized  to  achieve  the  high  gradients  needed.  The  RF 
cavity  is  a  metal  chamber  of  a  certain  shape  that  permits  particular  RF  modes  to  exist 
inside  of  it  by  matching  necessary  boundary  conditions.  However,  RF  cavities  come  with 
their  own  limiting  issues  that  must  be  addressed.  Ohmic  heating  in  the  cavity  walls 
caused  by  the  oscillating  electric  current  increases  the  amount  of  power  that  must  be 
pumped  into  the  system  to  accelerate  the  free  electron  bunch  and  could  require  that 
cooling  be  incorporated  into  the  structure  to  remove  this  heat.  Smaller  cavities  can  lead  to 
smaller  /  cheaper  accelerators,  but  as  the  cavity  size  shrinks,  the  resonant  frequency  of  the 
cavity  increases.  This  increased  frequency  increases  the  power  density  on  the  walls 
causing  more  heating  from  ohmic  losses.  To  counteract  this,  superconducting  RF  (SCRF) 
cavities  have  been  developed  to  decrease  the  ohmic  losses  to  nearly  zero. 

A  superconducting  material  has  a  resistance  that  goes  to  zero  for  DC  current  and 
very  close  to  zero  for  oscillating  currents  when  its  temperature  drops  below  a  critical 
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temperature.  These  SCRF  eavities  inerease  the  eomplexity  of  the  system  as  they  need  to 
be  eontained  inside  of  dual  layer  dewers  (cryomodules)  to  keep  them  below  this  eritieal 
temperature,  where  they  beeome  superconducting.  A  cryomodule  contains  two  major 
layers:  the  inner  vessel  contains  liquid  helium  in  direct  contact  with  the  outside  of  the 
cavity,  cooling  the  cavity  to  the  necessary  2  K  to  4  K  temperatures,  and  an  outer  vessel 
contains  a  liquid  nitrogen  bath  in  contact  with  the  liquid  helium  dewer  that  prevents 
excessive  boil  off  of  the  liquid  helium.  A  cryoplant  may  also  be  required  to  replace  the 
liquid  gases  that  have  evaporated,  which  adds  an  extra  component  to  the  system  and 
increases  complexity. 

For  the  Niowave  THz  FEL,  a  base  frequency  of  '-'350  MHz  is  being  used  in  a 
normal  conducting  copper  cavity  of  the  injector.  This  base  frequency  will  be  gated  to 
43.75  MHz  to  increase  the  charge  per  micropulse,  or  separated  group  of  electrons. 

2.  Linear  Accelerator 

After  the  electron  bunches  are  initially  accelerated  by  the  injector  they  are  not  up 
to  the  desired  energy  for  the  FEL  interaction;  therefore,  a  linac  is  formed  from  a  group  of 
RE  cavities  to  further  accelerate  the  beam.  The  number  and  size  of  the  accelerating 
cavities  depend  on  the  final  energy  and  the  RE  frequency  being  used;  the  cavities  can  also 
be  either  normal  conducting  or  superconducting.  Designs  for  linacs  generally  follow  a 
series  of  elliptical  or  cylindrical  cavities  placed  end-to-end.  An  example  can  be  seen  in 
Eigure  3. 
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Figure  3.  Cross  section  of  five  cell  superconducting  elliptical  cavity 

accelerator,  from  [25]. 

While  these  are  widely  used  shapes  for  both  normal  conducting  and  superconducting  RF 
accelerators,  they  are  not  the  only  ones  available.  Another  option  is  a  “spoke  cavity” 
configuration  to  accelerate  the  electron  beam.  This  geometry  contains  a  similar  cavity 
shape  with  a  “spoke”  in  the  middle.  This  spoke  has  an  axial  hole  in  it  that  acts  as  a  drift 
tube  where  the  electron  bunch  is  shielded  from  the  RF  fields  when  they  are  in  the  wrong 
phase  (as  seen  in  Figure  4). 


Figure  4.  Schematic  of  single  spoke  cavity,  from  [26]. 

Just  as  in  elliptical  cavity  linacs,  several  of  these  spoke  structures  can  be  placed  in  series 
to  provide  additional  acceleration.  One  main  benefit  of  this  geometry  is  that  for  the  same 
sustained  RF  frequency,  the  radius  of  a  spoke  cavity  is  approximately  one  half  that  of  an 
elliptical  cavity  [27],  so  a  lower  base  RF  value  can  be  used  without  increasing  the 
dimensions  of  the  cavity.  Lower  RF  values  are  desirable  since  they  decrease  ohmic 
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losses,  possibly  allowing  for  4  K  level  refrigeration  to  be  used  instead  of  a  larger,  more 
complex  2  K  cryoplant.  This  shape  also  allows  fewer  mechanical  vibrational  modes  to 
couple  to  the  EM  field  than  the  equivalent  elliptical  cavity  [28].  These  properties  make  a 
spoke  cavity  design  desirable  for  possible  shipboard  applications.  A  superconducting 
three-spoke  cavity  linac  is  chosen  for  the  Niowave  THz  FEE  to  produce  the  final  8  MeV 
electron  beam  energy. 

3.  Undulator 

An  undulator  is  needed  to  transfer  energy  from  the  relativistic  electron  beam  to 
the  light  pulse.  An  undulator  consists  of  two  rows  of  permanent  magnets,  each  with 
alternating  pole  faces,  on  either  side  of  the  beam  pipe.  This  creates  a  transverse  periodic 
magnetic  field  causes  the  electron  beam  to  wiggle.  The  Lorentz  force  causes  transverse 
oscillations  in  the  electron  bunch,  which  couples  the  electron  motion  to  the  optical  field 
generated  by  the  previous  electron  bunches,  and  allows  the  net  energy  transfer  to  occur. 
A  diagram  of  a  permanent  magnet  undulator  is  seen  in  Figure  5. 


FEL 


Electron  Beam 


Electron  Gun  + 
Linear  Accelerator 


Figure  5.  Permanent  magnetic  undulator  schematic  with  electron  beam,  from 

[29]. 


An  in-depth  analysis  of  this  interaction  will  be  discussed  in  a  later  chapter. 

4.  Electron  Beam  Control 

The  electron  beam  is  contained  in  vacuum  from  the  cathode  surface  to  the  beam 


dump  to  prevent  interaction  by  air  molecules.  The  electron  beam  needs  to  be  controlled 
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from  when  it  is  created  at  the  injector,  through  the  accelerating  sections  of  the  linac, 
directed  into  the  undulator,  and  then  finally  disposed  of  in  a  beam  dump.  This  steering  is 
done  with  static  magnetic  fields.  The  two  main  necessities  for  electron  beam  control  are 
the  ability  to  focus  the  electron  beam  and  to  bend  the  beam  to  stay  within  the  geometry  of 
the  beam  line.  A  set  of  two  or  three  quadrupoles  can  be  used  to  supply  focusing  in  both 
transverse  directions.  A  dipole  magnet  is  used  to  supply  a  magnetic  field  across  the 
electron  beam’s  path  to  alter  its  trajectory.  Design  of  these  magnets  will  be  discussed  in 
more  depth  in  a  later  chapter. 

5.  Optical  Cavity 

The  optical  cavity  surrounding  the  undulator  allows  the  optical  field  to  build  up 
over  multiple  passes  through  the  electron/light  interaction  region  in  the  undulator. 
However,  the  mirror  material  selected  can  decrease  the  ability  to  tune  the  PEL,  as  the 
mirrors  must  reflect  all  radiation  in  the  tunable  range.  For  the  Niowave  THz  FEE,  a  gold 
coated  mirror  has  been  chosen  due  to  its  high  reflectivity  in  the  THz  frequency  band. 
Further  design  considerations  will  be  discussed  in  a  later  chapter. 

One  common  method  uses  a  partially  transmissive  mirror  to  allow  a  small 
percentage  of  the  radiation  power  through.  This  is  often  a  preferable  method  as  it  does 
not  clip  or  distort  the  wavefront  of  the  optical  pulse,  but  it  can  only  be  utilized  for  certain 
wavelength  ranges  and  power  levels  due  to  its  dependency  on  the  optical  characteristics 
of  the  mirror  material  at  the  desired  wavelength.  Another  method  is  a  hole-outcoupling 
technique  where  the  mirror  has  a  small  hole  allowing  a  small  percentage  of  the  radiation 
to  escape.  The  choice  of  hole  outcoupling  for  the  Niowave  THz  optical  cavity  will  be 
discussed  in  a  later  chapter. 

B.  FEL  THEORY 

The  generation  of  coherent  radiation  stems  from  the  electron  bunches’  interaction 
with  the  magnetic  field  of  the  undulator  and  the  EM  fields  of  the  light. 
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1. 


Electron  Equations  of  Motion 


The  relativistic  Lorentz  force  between  an  electron  and  both  the  electromagnetic 
fields  of  the  undulator  and  the  optical  pulse  is  (in  cgs  units) 

p  =  = -e[E +  {p  xB)]  Eq.  2.1 

which  leads  to  an  acceleration 

-  —  (E  +  p  xB),  Eq.  2.2 

dt  mc^  ^  ^ 


where  m  is  the  electron  mass,  e  is  the  electron  charge  magnitude,  p  =  \ jc  is  the 
dimensionless  electron  velocity,  B  is  the  magnetic  field,  E  is  the  total  electric  field,  and 


is  the  relativistic  Eorentz  factor.  The  energy  equation  for  the  electrons 
is  the  fourth  component  of  the  relativistic  Eorentz  force  equation. 


dt 


y  =  -—P-e, 

me 


Eq.  2.3 


which  describes  the  change  in  energy  of  the  electron  subject  to  an  electric  field.  In  order 
for  energy  exchange  to  occur,  the  electron  must  have  a  component  of  the  velocity  p  in  the 
transverse  direction  to  couple  to  the  periodic  transverse  electric  field  of  the  light. 


2.  Electron  Motion  in  the  Undulator  Field 


The  relativistic  electron  experiences  the  largest  net  force  from  the  periodic 
magnetic  field  of  the  undulator.  Eor  simplicity  now,  it  is  assumed  that  the  transverse 
motion  of  the  electron  in  the  undulator  is  primarily  due  to  the  undulator  field,  and  so  that 
the  electron’s  energy  remains  approximately  constant  throughout  the  interaction.  A 
helical  undulator  has  a  magnetic  field  that  spirals  along  the  length  of  the  undulator 
according  to 

Bu  =  5o[cos(/Coz)  i  +  sin(/Coz) /]  ,  Eq.  2.4 

where  =  2n/Xo,  Xo  is  the  undulator  period,  and  Bo  is  the  undulator  field  strength. 
Substituting  this  equation  into  the  Eorentz  force  equation  (Eq.  2.2)  and  separating  p  into 
its  constituent  parts  leads  to 


15 


^  ^  sin(/Coz)  I  +  cos(koZ)j],  Eq.  2.5 

where  is  the  transverse  veloeity  of  the  eleetron.  Integrating  Eq.  2.5  into  perfeet  helical 

orbits  leads  to 

P±  =  --  [cos(/Coz)  I  +  sin(/Coz) /]  ,  Eq.  2.6 

gBA. 

where  K  =  is  the  dimensionless  undulator  parameter,  and  the  constants  of 

integration  are  set  to  zero. 

To  evaluate  a  small  energy  change  from  equation  (Eq.  2.3),  the  electric  field  of 
the  light  must  be  considered.  In  a  helical  undulator,  the  light  will  be  helically  polarized, 
with  an  electric  field  of  the  form  E  =  E[cos(xp)  i  —  sm(xp)]],  where  xp  =  kz  —  ayt  +  cp, 
k  =  2n / A,  A  is  the  optical  wavelength,  o)  =  2nf  =  2nclA  is  the  angular  frequency,  and 
(p  is  the  optical  phase.  Inserting  the  electric  field  into  Eq.  2.3  leads  to 

y  ^  Eq.  2.7 

Substituting  fix  and  fiy  from  Eq.  2.6  into  Eq.  2.7  leads  to 

.  _  —  (jjt  +  (o)  —  ^^cos(f  +  (o),  Eq.  2.8 

ymc  ^  ymc  v  ?  r  ^  ^  i 

where  ^  =  (/Cq  +  k)z  —  a>t  is  the  electron  phase,  which  describes  the  position  of  the 
electron  on  the  scale  of  an  optical  wavelength. 

3.  Resonance  Condition 

To  transfer  energy  from  an  electron  to  optical  field,  the  relation  y  <  0  must  be 
met.  Erom  Eq.  2.11,  y  oc  cos(^  +  cp),  loss  of  the  electron  energy  occurs  near  ^  +  cp  = 
(/Cq  +  k^p^ct  —  (ji)t  +  (p  ~  n.  To  ensure  that  the  electrons  persistently  lose  energy,  y  <  0 
and  ^  +  (p  ~  n  always,  the  electron  phase  must  not  change  significantly,  so  that 

(  =  (/Cq  +  k)PzC  —  (X)  ~  Q  .  Eq.  2.9 

Making  the  substitutions  that  co=kc,  ko=2n/Ao,  and  kF=2n/A  leads  to 

A  ^  Ao(i-fe)^  Eq.2.10 

HZ 
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which  describes  the  wavelength  of  light  that  has  the  maximum  energy  transfer  from  an 
electron  traveling  along  the  undulator  axis  with  speed  cf^z.  To  put  this  into  a  more  useful 
form,  a  few  simple  substitutions  ean  be  made.  The  definition  of  the  relativistie  Lorentz 

faetor,  =  1  —  |  —  1  —  combined  with  Eq.  2.6,  leads  to 

=  1-  [cos(/Coz)2  +  sm(ikoZy]-p/ .  Eq.  2.11 


Sinee  the  portion  inside  the  square  braekets  of  Eq.  2.1 1  is  identieally  equal  to  one. 


Pz=l- 


(I+K’2) 


Eq.  2.12 


2  2 

Expanding  the  square  root  as  a  Taylor  series,  and  for  relativistie  eleetrons  {\+K  )ly  «1, 
leads  to 


X 


Ao(i+/r^) 

2y2 


Eq.  2.13 


This  gives  the  approximate  wavelength  of  light  the  EEE  emit  at  for  a  partieular  undulator 
configuration  and  electron  beam  energy.  It  elearly  shows  that  EEEs  can  be  designed  and 
tuned  to  operate  at  any  wavelength,  depending  on  the  electron  beam  and  undulator 
parameters. 


4.  FEL  Pendulum  Equation 

It  is  essential  to  describe  the  EEE  interaction  by  looking  at  how  the  eleetron’s 
phase  evolves  as  the  eleetron  travels  through  the  undulator.  To  aid  in  this  approaeh,  a 
new  dimensionless  unit  of  time  is  introdueed,  t  =  - 1,  where  L=NXo  is  the  undulator 
length.  As  the  eleetron  traverses  the  undulator,  the  dimensionless  time  ehanges  from 

z=0->l.  This  leads  to  dr  —  -dt,  and  a  shorthand  notation  is  defined  as  i_)  =  d(_)ldT. 

Taking  the  derivative  of  the  eleetron  phase  with  respect  to  r  leads  to  the  electron  phase 
veloeity 

^  =  +  k)P,  -  k].  Eq.  2.14 

To  find  the  rate  of  ehange  of  the  eleetron  phase  veloeity,  v  is  differentiated  again  with 
respeet  to  r. 
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Eq.  2.15 


°  L  dv 

V  = - 

c  dt 


c 


(K  +  A 


c 


making  the  approximation  that  k»ko,  equivalent  to  assuming  y»\  near  resonanee. 
Taking  the  expansion  of  Eq.  2.12, 


(1+^A 

2  ^ 


Eq.  2.16 


and  then  taking  the  derivative  with  respeet  to  time, 


A 


-(i+zfA 


Eq.  2.17 


leads  to  an  expression  involving  the  ehange  in  eleetron  energy  (Eq.  2.8).  Recall  the 
definition  of  the  optical  wave  number,  which  at  resonance  (Eq.  2.13)  has  a  value  of 


k 


2y2 

Ao(l+^f2) 


Eq.  2.18 


Substituting  Eq.  2.18  and  Eq.  2.17  into  Eq.  2.15  leads  to 

V  =  — 4;/rN— , 
c  7 


Eq.  2.19 


where  N=L/Xo  is  the  number  of  undulator  periods.  Inserting  Eq.  2.8  into  Eq.  2.19  then 
gives 

V  =  |a|cos(i^  +  9),  Eq.  2.20 

where  |a|  =  ATiNeKL\E\/Y^Tnc^  is  the  dimensionless  optical  field  amplitude  ([A  is  the 
optical  electric  field  amplitude  in  cgs  units).  Eq.  2.20  is  the  same  form  of  the  differential 
equation  that  describes  the  evolution  of  a  simple  pendulum,  6  =  a  sin(0) .  This  allows  the 
EEE  physics  to  be  described  using  much  the  same  terminology  and  concepts  as  the 
simple  pendulum,  which  helps  gain  insight  into  the  interaction  between  the  electrons  and 
the  light.  Studying  the  phase  space  evolution  of  an  electron  during  the  EEE  interaction 
allows  an  intuitive  view  of  how  the  electron  energies  are  evolving  in  time  as  they  traverse 
the  undulator.  Eigure  6  shows  the  phase  space  diagram  for  a  simple  pendulum. 
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Pendulum  Phase  Space  Paths 


e 

Figure  6.  Pendulum  phase  spaee  plot. 

The  blue  line  denotes  the  “separatrix”  which  follows  the  phase  space  path  of  a 
pendulum  that  has  just  enough  energy  to  get  to  the  top,  6=n,  but  not  go  over.  All  orbits 
inside  of  this  line  are  “closed  orbits”  corresponding  to  pendula  swinging  back  and  forth 
with  amplitude  The  orbits  outside  of  the  separatrix  are  “open  orbits”  where  the 

pendula  have  enough  kinetic  energy  to  overcome  the  gravitational  potential  energy  and 
continually  swing  over  the  top,  moving  in  the  same  angular  direction.  The  pendulum’s 
orbit  is  determined  by  its  initial  conditions.  The  phase  space  trajectories  describe  the 
transfer  of  energy  between  the  potential  and  kinetic  energy  of  the  pendulum  and  gives 
insight  into  the  FEL  interaction. 

Translating  these  ideas  over  to  the  physics  of  an  FEL  interaction  involves  thinking 
about  the  electron  beam  as  a  whole.  As  was  stated,  the  electron  phase  measures  the 
position  on  the  scale  of  a  wavelength  of  light;  a  typical  electron  bunch  is  -1000  optical 
wavelengths  long,  and  has  -10^  electrons  in  each  optical  wavelength  X.  This  indicates  that 
the  electrons  are  initially  spread  randomly  over  each  optical  wavelength.  The  phase 
velocities  v  of  the  electrons  are  determined  by  the  electron’s  energies.  Eigure  7  shows 
snapshots  of  the  electron  beam  phase  space  as  it  evolves  through  the  undulator  from 
z=0-^  1  of  a  typical  EEL  on  resonance  with  low  gain.  The  blue  dots  represent  samples 
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electrons  while  the  separatrix  (red  line)  is  given  by 


of  2\a\ 


1/2 


^  =  2|a|  [1  +  sin(  +  )],  with  a  height 


A:  Electron  Phase  Space  Distribution:  x=0.0  B:  Electron  Phase  Space  Distribution: t=0. 2 


0  2  4  0  2  4 

C:  Electron  Phase  Space  Distribution:  x=0.4  D:  Electron  Phase  Space  Distribution: x=0. 6 


0  2  4  0  2  4 

E:  Electron  Phase  Space  Distribution: x=0. 8  A:  Electron  Phase  Space  Distribution: x=1 .0 


Figure  7.  The  phase  space  procession  of  the  electrons  during  the  FEL 

interaction  about  Vo=0. 


Figure  7  A  in  the  upper  left  shows  the  initial  phase  space  distribution  of  the 
electrons  at  t=0  with  a  phase  velocity  spread  about  Vo=0,  corresponding  to  resonance.  For 
the  sake  of  this  example  an  initial  field  amplitude  of  ao=25  was  used.  As  the  FEL 
interaction  progresses  (B,  C,  D,  E,  E)  the  electrons  follow  phase  space  trajectories  that  are 

described  by  the  EEL  pendulum  equation.  Electrons  at  an  initial  phase  (r>7r/2  lose  energy 
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(v  decreases),  while  electrons  at  an  initial  phase  C<7r/2  gain  energy  (v  increases).  By 
conservation  of  energy  the  only  place  for  the  energy  to  go  from  the  electrons  is  into  the 
optical  field.  In  this  example,  about  half  the  electrons  gain  energy  and  half  lose  energy,  so 
very  little  net  energy  is  transferred  into  the  optical  field.  Figure  8  shows  how  to  make  the 
FEL  “work”  by  operating  off  resonance. 


A:  Electron  Phase  Space  Distribution: 'c=0.0 


C:  Electron  Phase  Space  Distribution:  t:=0.4 


E:  Electron  Phase  Space  Distribution:  ^=0.8 


c 


B:  Electron  Phase  Space  Distribution:  ^=0. 2 

15| - ' - ' - ^ - ' - ' - ' - 


-10  12  3  4 

c 


D:  Electron  Phase  Space  Distribution:  t=0. 6 


A:  Electron  Phase  Space  Distribution:  t^I  .0 


^ _ 

Figure  8.  The  phase  space  procession  of  the  electrons  during  the  FEL 

interaction  about  Vo=6. 

As  in  Eigure  7,  Eigure  8  shows  the  progression  of  the  phase  space  trajectories  for 
an  electron  beam  as  it  undergoes  the  EEL  process.  In  this  case,  the  initial  electron  phase 
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velocity  is  centered  off  resonance  around  Vo=6.  This  leads  to  the  majority  of  electrons 
who  start  inside  the  separatrix  to  lose  energy  as  z=0^1.  This  net  energy  transfer  leads  to 
an  increase  in  the  optical  field. 

5.  Optical  Field  Evolution 

For  the  FEL  process  to  be  understood  self-consistently,  it  is  necessary  to  study  the 
optical  field  evolution. 

a.  Spontaneous  Emission 

Light  is  initially  created  by  spontaneous  emission  from  the  transverse  acceleration 
of  the  relativistic  electron  beam.  This  radiation  is  projected  in  a  forward  cone  with  an 
angle  width  of  ^Hy.  The  spontaneous  emission  will  have  a  spectral  width  of  -'1/A^ 
centered  on  the  resonant  wavelength  of  the  FEL. 

b.  Weak  Field  Gain 

Eor  the  EEL  process  to  work,  net  energy  must  be  transferred  from  the  electron 
beam  to  the  optical  pulse.  Erom  the  average  electron  phase  velocity  change  of  the 
electron  beam,  (Av),  one  can  determine  the  average  loss  of  energy  by  the  bunch, 
(Aymc^),  from  Eq.  2.22,  and  thus  the  energy  gained  by  the  optical  field.  As  the  EEL 
turns  on,  the  optical  fields  produced  by  spontaneous  emission  are  very  small,  |a|«7i.  In 
weak  fields  with  low  gain,  \a\~ao.  Elsing  perturbation  theory,  one  can  solve  the  pendulum 
equation  in  powers  of  Uq.  To  0**^  order  in  Uo,  the  solution  is  trivial:  with  no  optical  field 
there  is  no  energy  change.  To  order  in  Uo,  when  the  solution  is  averaged  over  the 
whole  electron  beam  spread  over  an  optical  period  the  resultant  energy  change  is,  again, 
identically  zero.  One  must  go  out  to  2"‘*  order  in  Oo  before  finding  a  non-zero  solution 
[30], 

(Av)  =  -^^4"  [cos(Vot)  —  1+  VoTsin(VoT)/2]  .  Eq.  2.21 
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To  find  the  total  energy  ehange  of  the  eleetron  beam,  AU,  one  must  substitute  Eq.  2.21 
into  Eq.  2.19  and  multiply  by  the  total  number  of  eleetrons  in  the  beam,  ne=pVe,  where  p 
is  the  electron  beam  density  and  Ve  is  the  volume  of  the  electron  bunch 

The  EEL  gain  is  defined  as  the  change  in  energy  of  the  optical  pulse  with  respect 
to  the  initial  energy  in  the  optical  pulse.  The  energy  in  the  optical  field  is  given  by 
Uopt=E  VopJ^n,  where  E  is  the  electric  field  amplitude  and  Vopt  is  the  volume  of  the 
optical  pulse  (cgs  units).  Taking  all  of  this  into  account  leads  to  an  equation  for  the  weak 
field  gain,  G, 

^  _  (ymG) 

N 

and  using  Eq.  2.21, 

G  =  ^  [2  —  2  cos(Vot)  —  VqT  sin(VoT)],  Eq.  2.23 

'^o 

where  J  =  Sn^Ne^K^L^p/Y^nT-C^  is  the  dimensionless  current  density,  and  F=Ve/Vopt 
=Ae/Aopt  is  the  filling  factor,  which  denotes  how  well  the  electron  beam  with  area  Ae  and 
the  optical  mode  with  area  Aopt  overlap.  This  solution  is  valid  for  weak  fields  (a«7i)  and 
for  low  gain  (/«1).  Round  trip  gain  must  be  higher  than  the  round  trip  losses  in  the 
optical  cavity  for  the  EEL  process  to  work.  The  amount  of  gain  per  pass  is  dependent  on 
many  factors,  including  the  strength  of  the  optical  field,  the  initial  electron  phase 
velocity,  and  the  current  density.  Eigure  9  shows  the  weak  field  gain  as  a  function  of  Vq. 


(^v)pVe 


E^V, 


o  ^ovt J 


Eq.  2.22 
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Weak  Field  Gain 


Figure  9.  Weak  field  gain  as  a  funetion  of  Vq. 


The  gain  eurve  is  approximately  anti-symmetric  about  Vo=0  with  the  peak  gain  at  an 
initial  electron  phase  velocity  Vo~2.6.  If  Vo~-2.6,  then  the  gain  is  negative  and  at  a 
minimum,  meaning  that  the  electron  bunch  will  take  the  maximum  amount  of  energy 
from  the  optical  pulse.  At  Vo=0  the  gain  is  zero,  meaning  that  overall  the  electron  bunch 
neither  gains  nor  loses  energy  to  the  optical  pulse,  which  recreates  the  same  result  that  is 
achieved  from  the  analysis  of  the  phase  space  evolution  shown  in  Figure  7.  The  gain 
spectrum  is  approximately  anti-symmetric  with  a  gain  bandwidth  of  zlvo~47i.  From  Eq. 
2.19  the  gain  bandwidth  in  energy  is  approximately  UN. 

c.  Coherence  Development 

Like  other  lasers,  FELs  can  have  a  relatively  narrow  spectrum.  The  light  created 
from  spontaneous  emission  is  relatively  broadband,  corresponding  to  a  wide  range  of 
initial  electron  phase  velocities  centered  on  resonance  (vo=0).  As  the  EEL  process 
continues,  energy  is  preferentially  transferred  from  the  electron  beam  to  the  optical  field 
at  wavelengths  near  the  peak  of  the  gain  curve  (vo=2.6).  As  the  optical  pulse  bounces 
back  and  forth  in  the  cavity,  gaining  energy  from  each  subsequent  electron  bunch,  the 
optical  spectrum  becomes  narrower  and  narrower  (see  figure  10),  eventually  becoming 
transform  limited. 
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Figure  10.  Coherence  development  over  2500  passes,  from  [31]. 


Figure  10  depicts  the  power  spectrum  of  the  optical  pulse  in  an  FEL  over  2500  passes.  As 
the  number  of  passes  (n)  increases  the  width  of  the  spectrum  narrows  as  the  gain 
spectrum  preferentially  amplifies  wavelengths  near  peak  gain. 


d.  Wave  Equation 


While  the  pendulum  equation  describes  how  the  electrons  evolve  during  the  FEL 
process,  the  wave  equation  is  needed  to  describe  how  the  light  field  evolves  through  the 
undulator.  To  simplify  the  wave  equation,  the  “slowly-varying  envelope  approximation” 
(SVEA)  is  used,  which  states  that  the  phase  and  amplitude  of  light  do  not  change  much 
over  one  oscillation  of  the  electromagnetic  wave.  This  is  a  reasonable  approximation  for 
most  lasers  and  is  commonly  used.  With  this  framework,  the  wave  equation  can  be 
written  in  terms  of  the  vector  potential,  A(x,t),  and  the  transverse  current  density. 


Eq.  2.24 


A  solution  of  the  form 


A{x,t) 

k. 


Eq.  2.25 
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can  be  assumed,  whieh  corresponds  to  a  helieally-polarized  plane  wave  produeed  by  a 
helical  undulator.  Breaking  this  equation  into  its  spatial  and  temporal  eomponents  leads 
to  a  simpler  solution.  Assuming  no  transverse  dependenee  on  the  spatial  eomponent  (a 
plane  wave  where  E=E{z,t))  leaves  only  a  longitudinal  component  (d^/dz^)  in  the 
Laplacian.  This  is  a  good  approximation  for  a  beam  with  little  diffraetion  during 
evolution  through  the  undulator.  This  leads  to 


d^A 

dz^ 


k 


i-ii+j). 


Eq.  2.26 


where  the  second  derivative  of  E  is  ignored  due  to  the  SVEA.  The  temporal  term  ean  be 


written  in  mueh  the  same  manner: 


_  »  (-2.6,-  -  6.2£ j— +;). 


Eq.  2.27 


Plugging  Eq.  2.27  and  Eq.  2.26  into  Eq.  2.24  and  multiplying  by  the  eomplex  eonjugate 


of  the  unit  veetor  leads  to 


Eq.  2.28 


A  further  development  is  to  follow  a  slice  of  the  eleetric  field  as  the  light  travels  through 
the  undulator  at  speed  c.  This  has  the  benefit  of  removing  the  z  dependence  on  the  left 
hand  size  of  Eq.  2.28.  The  current  density  is  a  sum  of  individual  electron  trajectories  and 
can  be  written  out  as  /x  =  ~  where  Ej  is  the  position  of  the  z* 

electron.  The  electrons  follow  the  helical  trajectories  derived  in  Eq.  2.9,  so  the  eurrent 


ean  be  written  out  as 


/x  =  ecK  Re[i(-ii  +  — 5®(x-rj)] 

Yi 


Eq.  2.29 


^  g-ikoz 

/x  ~  ecKpRe[i(-ii  +  ;))  i—^)]  , 


Eq.  2.30 


where  the  average  over  the  eleetrons  in  the  beam,  (...),  has  been  taken,  and  p  is  the 
electron  beam  density.  Combining  Eq.  2.30  with  Eq.  2.28  then  leads  to  the  equation 


dE  g-il(ko+k)z-o)t] 

-  =  -2necKp  ( - - >  . 


Eq.  2.31 


This  describes  how  the  optieal  field  evolves  through  the  undulator  during  the  EEE 


proeess.  Transforming  this  into  dimensionless  units  and  assuming  that  all  the  electrons  in 
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the  bunch  have  approximately  the  same  energy,  y~yo,  leads  to  the  dimensionless  PEL 
wave  equation, 

a  =  -jF(^e~‘^),  Eq.  2.32 

where  a  and  j  are  the  dimensionless  optical  field  and  the  dimensionless  current  density  as 
previously  defined.  The  bracketed  term  is  known  as  the  “bunching”  term  and  denotes 
how  much  the  electrons  have  bunched  on  the  order  of  a  wavelength  of  light.  If  the 
electrons  are  equally  spaced  out  over  the  wavelength  of  light  then 

=  0  ,  Eq.  2.33 

so  there  is  no  change  to  the  optical  field.  It  is  only  after  some  transit  time  through  the 
undulator  where  some  electrons  have  lost  energy  and  some  have  gained  energy  that  they 
start  to  bunch  on  the  scale  of  the  wavelength  of  light.  It  is  this  bunching  at  a  certain 
optical  phase  which  leads  to  energy  transfer  to  the  optical  field. 

e.  Short  Pulse  Effects 

Until  now  the  electron  bunch  length  has  not  played  a  role.  In  real  life,  the  electron 
bunch  is  of  finite  length  (often  ~ps),  which  must  be  taken  into  consideration. 

(1)  Optical  Eethargy 

As  was  previously  discussed,  the  electrons  must  be  bunched  on  the  scale  of  a 
wavelength  of  light  in  order  to  have  net  energy  transfer  between  the  electron  bunch  and 
the  optical  field.  It  takes  some  time  for  this  bunching  to  develop  as  the  electrons  interact 
with  the  optical  field  and  the  undulator  magnetic  field,  so  energy  transfer  does  not  happen 
at  the  very  beginning  of  the  undulator  but  rather  toward  the  end.  A  depiction  of  this  can 
be  seen  in  Eigure  1 1 . 
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Figure  1 1 .  Depiction  of  optical  lethargy  of  the  light  pulse  transiting  the 

undulator. 

In  Figure  11,  the  electron  bunch  (red)  and  the  optical  pulse  (blue)  are  co¬ 
propagating  through  the  undulator  (green  arrows)  from  left  to  right.  At  the  beginning  of 
the  undulator  the  electron  bunch  is  slightly  ahead  of  the  optical  pulse.  Although  the 
electrons  are  highly  relativistic,  they  still  are  not  traveling  at  the  speed  of  light,  so  the 
optical  pulse  overtakes  it  during  transit.  Halfway  through  the  undulator  the  electrons  are 
becoming  bunched  so  more  energy  is  transferred  to  the  optical  field  and  its  amplitude 
begins  to  grow.  By  the  end  of  the  undulator  (right),  the  electrons  are  well  bunched  and 
transfer  the  most  energy  to  the  optical  pulse.  This  has  the  effect  of  amplifying  the  rear  of 
the  optical  pulse  more  than  the  front,  causing  the  centroid  of  the  optical  pulse  to  travel 
less  than  the  speed  of  light.  This  effect  must  be  taken  into  consideration  when  designing 
the  FEL  oscillator  cavity. 

(2)  Desynchronism 

The  cavity  length  of  the  oscillator  must  be  properly  chosen  so  that  the  optical 
pulse  is  synchronized  with  each  subsequent  electron  bunch.  This  distance  is  based  on  the 
repetition  rate  of  the  electron  bunches  and  the  speed  of  light,  but  as  was  previously 
discussed  the  centroid  of  the  optical  pulse  is  less  than  c.  In  order  to  counteract  this  effect, 
the  length  of  the  cavity  must  be  shortened  slightly  (on  the  order  of  microns)  to  allow 
proper  overlap  between  the  electron  bunches  and  optical  pulse  over  many  passes.  If  this 
were  not  taken  into  account  then  the  tail  of  the  pulse  will  continue  to  grow  more  than  the 
front  of  the  pulse  and  eventually  would  evolve  away  from  the  electron  beam  and  the  FEL 
would  not  work.  Optimal  values  to  shorten  the  oscillator  cavity  to  account  for 
desynchronism  are  typically  between  zero  and  0.25  NX. 
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/.  Electron  Beam  Quality 


An  pel’s  performance  is  affected  by  the  quality  of  the  electron  beam.  A  key 
metric  of  electron  beam  quality  is  the  emittance,  which  is  defined  as  the  amount  of 
volume  the  electron  bunch  takes  up  in  phase-space.  There  are  two  types  of  emittance; 
transverse  emittance  denotes  how  the  electron  beam’s  transverse  size  and  angular  spread 
are  correlated  and  is  given  by  =  OyOg,  where  Or  is  the  rms  radius  of  the  beam  and  ae  is 
the  rms  spread  in  angles.  Longitudinal  emittance  is  related  to  the  energy  spread  {AE)  and 
the  pulse  length  of  the  electron  bunch  {Atrms)  and  is  given  by  £i  =  AE  ■  The  ideal 

electron  beam  for  the  LEL  interaction  would  be  monoenergetic  and  perfectly  collimated 
so  that  all  the  electrons  could  match  the  resonance  condition,  but  this  is  not  feasible. 
There  is  a  trade  off  with  electron  beam  generation  between  beam  quality  and  the  amount 
of  current  that  can  be  achieved  from  an  electron  accelerator  that  must  be  balanced 
depending  on  the  needs  of  the  LEL. 

Poor  electron  beam  quality  degrades  the  LEL  performance  by  decreasing  the 
amount  of  bunching  that  can  occur.  As  energy  is  directly  related  to  the  electron  phase 
velocity,  an  energy  spread  can  be  visualized  as  a  spread  in  Vo  (see  Ligure  12). 
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Figure  12.  Comparison  of  energy  spread  of  initial  electron  bunches,  strong 

fields. 


Figure  12  depicts  simulations  of  an  electron  bunch  with  an  initial  optical  field  of 
Uo  =  25,  7=30  with  two  different  energy  spreads.  The  initial  electron  phase-space 
diagrams  of  a  good  electron  bunch  with  a  fairly  minimal  energy  spread  of  0.5  percent 
(left)  as  well  as  one  with  a  larger  two  percent  energy  spread  (right)  can  be  seen  on  the  top 
row.  The  lower  two  plots  show  the  phase  space  diagrams  at  the  end  of  the  undulator  after 
the  FEL  interaction.  The  quality  of  the  bunching  is  lower  with  higher  energy  spread, 
which  decreases  the  gain  of  the  FEL.  The  effect  can  be  greater  in  the  weak  field  gain 
regime  where  electrons  with  a  large  phase  velocity  distribution  are  not  well  trapped  by 
the  weak  fields.  An  example  of  this  can  be  seen  in  Eigure  13. 
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Figure  13.  Comparison  of  energy  spread  of  initial  electron  bunches,  weak 

fields. 

In  the  weak  field  regime  with  a  small  initial  field  (ao=3  in  Figure  13),  the  large 
two  percent  energy  spread  does  not  allow  bunching  to  occur  and  reduces  the  weak  field 
gain.  This  is  a  key  reason  to  reduce  the  energy  spread  of  the  electron  beam. 

The  transverse  emittance  has  the  same  effect  of  increasing  the  spread  in  initial 
phase  velocities.  A  spread  in  angles  is  equivalent  to  a  spread  in  longitudinal  velocities. 
Since  is  proportional  to  Pz,  this  induces  a  spread  in  phase  velocities  in  addition  to  the 
initial  energy  spread,  which  again  decreases  the  bunching  and  hence  the  output  power. 
Likewise,  a  larger  beam  radius  will  lead  to  stronger  betatron  focusing,  thus  again 
increasing  the  phase  velocity  spread  and  decreasing  bunching  and  extraction. 
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III.  FEL  SIMULATIONS 


A.  NOMINAL  FEL  PARAMETERS 

The  nominal  parameters  simulated  for  the  Niowave  FEL  are  listed  in  Tables  2-4. 
The  electron  beam  energy  is  8  MeV,  the  bunch  charge  is  0.2  nC,  and  the  beam  radius  at 
waist  is  0.45  mm  (rms).  The  assumed  pulse  shape  was  parabolic  with  a  5  ps  full  width 
with  an  8.75  mA  average  current.  It  has  a  normalized  emittance  of  28  mm-mrad,  an 
energy  spread  of  1.66%,  and  a  repetition  rate  of  43.75  MHz. 


Beam  kinetic  energy 

8  MeV 

Bunch  charge 

0.2  nC 

Beam  radius  at  waist,  rms 

0.45  mm 

Longitudinal  pulse  shape 

Parabolic 

Pulse  duration,  full  width 

5  ps 

Average  current 

8.75  mA 

Normalized  emittance,  rms 

28  mm-mrad 

Energy  spread  (Ay/y) 

1.66% 

Repetition  rate 

43.75  MHz 

Table  2.  Electron  beam  parameters  being  simulated. 


To  increase  the  peak  current  to  60  A,  the  injector  will  be  gated  to  the  eighth  sub¬ 
harmonic  of  the  base  350  MHz  frequency  or  43.75  MHz  [32]. 


Undulator  period  (/lo) 

3.3  cm 

Number  of  periods  (N) 

10 

Undulator  parameter  (K),  rms 

0.72 

Undulator  length  (L) 

33  cm 

Table  3.  Elndulator  parameters  being  simulated. 
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The  undulator  has  A^=10  periods  of  Xo=2).3  em  leading  to  a  33  em  length  and  a  variable 
gap  that  ean  be  adjusted  to  supply  a  magnetic  field  corresponding  to  an  undulator 
parameter  range  of  0.72<if<1.42.  This  has  a  direct  effect  on  the  resonant  radiation 
wavelength  and  thus  the  optical  cavity  parameters.  The  resonant  wavelength  and  the 
optical  beam  waist  radius  determine  the  geometry  for  the  optical  cavity.  The  mirrors’ 
radius  of  curvature  ('-'1.7  m),  the  cavity  length  (-'3.42  m),  the  nominal  radiation 
wavelength  (93  pm)  determine  the  Rayleigh  length  Zo=ll  cm  for  the  cavity  leading  to  a 
mode  waist  radius  of  lFo=l-82  mm.  The  mirror  radius  is  eight  centimeters.  The  round  trip 
loss  was  assumed  to  be  6%  due  to  output  coupling  and  mirror  losses.  These  parameters 
are  listed  in  Table  4. 


Cavity  length  (S) 

3.428  m 

Rayleigh  length  (zo) 

11  cm 

Resonant  wavelength  (k) 

93  pm 

Mode  waist  radius  (1/e) 

1.82  mm 

Mirror  radius  of  curvature 

-1.7  m 

Mirror  radius 

8  cm 

Round  trip  losses 

6% 

Table  4.  Optical  cavity  parameters  being  simulated. 


B.  CODE  DESCRIPTIONS 

Computer  codes  have  been  developed  that  model  the  PEL  interaction  and  are 
based  on  the  theory  first  described  in  the  1970s  [33].  Three  different  FEE  codes  are 
discussed  and  used  in  this  dissertation:  the  first  is  GINGER,  a  widely  used  code  created 
by  scientists  at  Eawrence  Berkeley  National  Eaboratory  [34];  the  second  is  4DPULSEN 
developed  by  the  Naval  Postgraduate  School;  the  last  is  a  code  was  written  by  the  author 
of  this  dissertation. 
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1. 


GINGER 


GINGER  is  a  transverse  multidimensional  (ineluding  r-t  components  in  the  fields 
and  x-y-t  in  the  electrons)  and  polychromatic  PEL  simulation.  For  this  study  the  code  was 
run  in  a  monochromatic  mode,  which  ignores  longitudinal  pulse  shape  effects  but  does 
take  into  account  optical  beam  diffraction.  One  benefit  of  GINGER  is  its  ability  to  input 
electron  beam  parameters  from  external  particle  tracking  codes  relatively  easily. 

2.  4DPULSEN 

Researchers  at  the  Naval  Postgraduate  School  have  created  a  suite  of  FEE  codes 
that  model  the  interactions  in  a  number  of  ways.  One  of  the  newest  is  a  4D  pulse  code 
that  models  the  electron  beam  and  the  optical  field  in  x,y,z,  and  t,  and  includes  pulse 
slippage,  desynchronism,  and  diffraction  effects.  The  code  can  model  an  FEE  oscillator 
and  is  used  to  study  stability  criteria  by  shifting  and  tilting  the  mirrors.  The  code  is 
parallelized  and  runs  on  a  cluster  computer;  with  ~100  nodes  the  simulation  takes  ~3 
seconds  per  pass,  and  may  run  over  thousands  of  passes. 

3.  ID  EEL  Code 

An  FEE  code  was  written  by  the  author  of  this  dissertation  to  calculate  gain  and 
extraction  (percentage  of  electron  beam  power  transferred  to  the  optical  pulse).  This  is  a 
ID  code  that  tracks  the  electrons’  phase  space  evolution  as  they  traverse  the  undulator 
using  the  FEE  pendulum  and  wave  equations  discussed  in  Chapter  2.  It  utilizes  the  Euler- 
Cromer  method  to  iteratively  step  through  the  coupled  partial  differential  equations.  In 
each  time  step  At,  the  change  in  each  electron’s  phase  velocity  is  calculated  using 
Av  =  (Uj.  cos(^)  —  ajSin(^))At,  where  Ur  and  a,  are  the  real  and  imaginary  parts  of  the 
dimensionless  optical  field.  Next  the  new  electron  phase  velocity,  v(t  +  At)  =  v(t)  + 
Av,  is  used  to  update  the  electron  phase,  ((t  +  At)  =  (^(t)  +  v(t  +  At) At.  This  is  done 
for  each  individual  electron.  At  the  end  of  the  time  step,  the  changes  in  the  optical  field  is 
calculated  using  the  expanded  version  of  Eq.  2.34,  leading  to 

AUj.  =  —jF{cos(Q)  ■  At ,  Eq.  3.1 
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and 


Aui  =  jF {sm(Q)  ■  At, 


Eq.  3.2 


where  (...)  represents  an  average  over  the  sample  eleetrons  and  j  is  the  dimensionless 
eurrent  density  and  F  is  the  filling  factor.  This  new  optical  field  is  then  used  to  calculate 
the  electron  phase  velocity  change  Av  of  the  next  time  step.  In  updating  the  electron 
phases,  periodic  boundary  conditions  are  used  to  maintain  -%I2<C<2)%I2  for  each  electron. 
The  initial  optical  power,  +  aj^^,  is  calculated  at  the  beginning  of  the 

undulator  from  the  initial  components  Uro  and  Uio.  The  real  and  imaginary  components  at 
the  end  of  the  undulator  {arf  and  ay  respectively)  are  used  to  calculate  the  final  power, 
.  The  gain  is  calculated  from  the  fractional  change  in  the  optical  power 
over  the  initial  optical  power. 


G  = 


Eq.  3.3 


Key  inputs  into  the  code  include  the  initial  dimensionless  field  components  a^  and  a/o, 
the  dimensionless  current  density  j,  the  filling  factor  F,  and  the  energy  spread  of  the 
electron  beam.  This  code  is  useful  for  its  ability  to  show  trends  in  the  performance  of  the 
EEE  by  being  able  to  do  a  large  number  of  simulations  quickly  and  on  a  machine  with 
minimal  processing  power. 


C.  GAIN  SPECTRUM 


The  weak  field  gain  (G)  of  the  EEE  describes  how  the  optical  pulse  gains  energy 
over  each  pass  before  the  onset  of  saturation.  The  gain  is  dependent  on  the  initial  electron 
phase  velocity  (amongst  other  things).  Eigure  14  compares  the  ID  simulation  results  to 
the  analytic  solution  (Eq.  2.26). 
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Gain  Spectrum 


Figure  14.  ID  simulation  results  (blue)  and  analytic  solution  (red)  of  gain 
spectrum  of  Niowave  THz  FEL  for  0%  energy  spread. 

The  analytic  solution  does  not  take  into  account  energy  spread,  so  the  spread  was  set  to 
zero  for  this  simulation;  all  other  parameters  are  the  nominal  parameters  in  Tables  2-4 
corresponding  to  a  dimensionless  current  density  j~3.  The  ID  simulation  (blue)  obtains  a 
maximum  G~46%  at  Vo  ~2.5,  similar  to  the  analytic  solution  (red).  The  self-consistent 
update  of  the  radiation  field  in  the  simulation  (blue)  causes  the  asymmetry  of  slightly 
smaller  loss  compared  to  the  analytic  solution  where  the  field  is  held  constant.  Flowever, 
care  must  be  taken  when  analyzing  these  results.  Diffraction  and  longitudinal  pulse 
effects  of  the  FEL  will  cause  a  decrease  in  the  gain.  The  slippage  length  (M~l  mm)  is  on 
the  order  of  the  pulse  length  (5  ps->~l  mm),  which  will  decrease  the  overlap  of  the 
electron  beam  with  the  optical  pulse. 

The  ID  code  can  also  be  used  for  non-monoenergetic  beams.  An  energy  spread 
corresponds  to  a  spread  in  phase  velocities,  which  will  cause  different  electrons  to  be 
spread  into  different  portions  of  the  gain  spectrum.  A  depiction  of  this  can  be  seen  in 
Eigure  15. 
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The  graphs  in  Figure  15  show  the  same  nominal  gain  speetrum  (blue)  that  was 
simulated  in  Figure  14  with  an  overlay  of  a  nominal  energy  spread  in  the  eleetron  beam 
(red).  The  left  graph  shows  a  small  energy  spread  where  approximately  all  eleetrons  still 
eontribute  to  positive  gain  eausing  the  peak  gain  to  still  oecur  at  Vo~2.6.  In  the  right  graph 
of  Figure  15,  the  energy  spread  is  mueh  larger  resulting  in  some  eleetrons  experieneing 
reduced  or  even  negative  gain.  Some  of  the  electrons  on  the  left  side  of  the  distribution 
experience  negative  gain,  taking  energy  out  of  the  optical  pulse.  The  middle  portion  still 
experience  positive  gain  while  those  to  the  right  of  the  peak  have  reduced  gain.  The  net 
effect  causes  the  value  of  Vo  giving  the  best  gain  to  shift  to  the  right.  Applying  this  code 
to  the  realistic  electron  beam  of  the  Niowave  TFIz  FEL  with  an  energy  spread  of 
Jy/y~l  .66%  leads  to  the  resultant  gain  spectrum  in  Figure  16. 
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Figure  16.  Gain  spectrum  of  Niowave  TFIz  FEL:  Ay/7=1.66%. 

The  gain  spectrum  has  the  same  approximate  anti-symmetry  as  the  nominal  mono- 
energetic  beam  case  but  some  new  key  features  are  that  the  peak  gain  is  reduced  from 
-46%  to  -24%,  the  position  of  the  peak  is  shifted  slightly,  from  Vo~2.6  to  Vo~3  for  the 
reasons  discussed  above,  and  the  side  peaks  of  local  maxima  have  disappeared.  Recall 
that  a  shift  in  Vo  corresponds  to  a  shift  in  wavelength,  so  this  modified  gain  spectrum 
alters  the  wavelength  that  is  preferentially  amplified  on  each  pass  through  the  undulator. 

Using  the  parameters  in  Tables  2-A,  GINGER  can  be  used  to  take  into  account 
energy  spread  and  diffraction  effects  of  the  optical  pulse.  The  resulting  gain  spectrum  is 
shown  Figure  17. 


39 


Figure  17.  Scan  of  weak  field  gain  vs.  wavelength  using  GINGER  FEL  code  for 

the  Niowave  FEL  parameters. 

The  shape  shows  good  agreement  with  Eigure  16  including  the  approximate  anti¬ 
symmetric  shape  of  the  gain  spectrum  and  the  flattening  of  the  side  lobes  due  to  the 
energy  spread,  and  approximately  the  same  peak  gain  of  G~22%. 

D.  WEAK  FIELD  GAIN  VS.  ENERGY  SPREAD 

The  gain  was  studied  as  a  function  of  energy  spread  in  the  electron  beam  using 
the  three  codes.  Simulations  were  done  varying  the  energy  spread  from  zero  percent,  the 
ideal  case,  up  to  3.25  percent,  an  extreme  case.  Eor  all  simulations,  a  Gaussian  energy 
spread  shape  was  assumed.  The  results  from  the  ID  code  can  be  seen  in  Eigure  18. 
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Figure  18.  Gain  degradation  from  energy  spread  using  ID  code. 

With  no  energy  spread  the  ID  simulations  show  the  FEL  would  have  G~  46  and 
in  the  worst-case  scenario,  a  3.25  percent  energy  spread  leads  to  G~10%.  For  the 
Niowave  THz  FEL,  an  energy  spread  of  Ay/y~\.7%  is  expected  from  the  linac  from 
preliminary  measurements  and  simulations  corresponding  to  a  G~25%  in  Eigure  18.  The 
green  line  denotes  the  expected  six  percent  round  trip  losses  of  the  cavity.  While  this  is  a 
relatively  high  gain  and  allows  for  a  large  amount  of  loss,  this  simple  ID  code  does  not 
take  into  account  diffraction  or  longitudinal  pulse  effects. 

To  analyze  the  effects  of  diffraction  and  short  pulses  on  the  expected 
performance,  GINGER  and  4DPULSEN  can  be  used  to  scan  through  the  energy  spread. 
These  results  are  seen  in  Eigure  19. 
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Figure  19.  GINGER  (red)  and  4 DPULSEN  (bluo)  eode  results  of  gain  against 
initial  eleetron  energy  spread  for  the  Niowave  THz  FEL. 

The  results  for  GINGER  (red  line),  whieh  takes  into  aeeount  diffraetion,  show  a 
redueed  gain  of  G~32%  for  the  ideal  ease  and  a  value  of  G~22%  at  the  expeeted  energy 
spread.  As  was  stated,  in  the  monoehromatie  mode  that  was  used,  GINGER  does  not  take 
into  aeeount  the  longitudinal  pulse  profile  of  the  electron  beam  and  is  therefore 
optimistic.  Short  pulse  effects  will  cause  the  FEL  gain  to  decrease  due  to  poor  overlap 
between  the  radiation  pulse  and  electron  pulse  during  the  interaction,  as  can  be  seen  by 
the  reduced  gain  the  4DPULSEN  code  predicts  (blue).  This  code  predicts  only  a  G=22% 
for  the  monoenergetic  beam  case  and  G~13%  for  the  expected  energy  spread.  Eor  each 
data  point  of  the  4D  code  the  amount  of  desynchronism  was  altered  to  maximize  the  gain. 
A  few  results  of  the  desynchronism  scan  can  be  seen  in  Eigure  20. 
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Figure  20.  4DPULSEN results  showing  gain  vs.  dimensionless  desynchronism 
(d)  for  an  energy  spread  of  0%  (blue),  1.2%,  (red),  and  2.4%  (black). 

The  optimum  desynchronism  decreases  as  the  energy  spread  increases.  This  shift 
of  the  optimum  desynchronism  d  corresponds  to  decrease  in  the  cavity  length  of  ~13  pm. 
The  4DPULSEN  co&e  results  of  only  G~13%  for  the  expected  energy  spread  in  Figure  19 
will  be  the  assumed  value  for  the  experiment  will  be  built  around  to  ensure  that  the  round 
trip  losses  are  less  than  the  gain. 

E.  VARYING  UNDULATOR  GAP  HEIGHT 

Adjusting  the  transverse  gap  height  of  the  undulator  changes  the  peak  magnetic 
field,  which  alters  the  FEL  resonant  wavelength  (all  other  things  being  equal).  This 
change  to  longer  wavelengths  causes  the  optical  mode  to  expand  at  a  different  rate  due  to 
diffraction.  This  diffraction  needs  to  be  taken  into  consideration  when  designing  the 
optical  cavity  and  beam  line,  as  extensive  clipping  will  cause  the  round  trip  losses  to  be 
greater  than  the  gain.  As  the  undulator  gap  width  is  decreased,  the  peak  magnetic  field 
gets  stronger,  increasing  K,  causing  the  resonance  wavelength  to  increase.  In  terms  of  the 
cavity  design  this  is  detrimental  for  two  reasons:  as  the  gap  width  is  decreased,  the 
clipping  on  the  beam  pipe  inside  the  undulator  increases,  and  a  narrower  gap  corresponds 
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to  longer  wavelength,  which  leads  to  more  diffraction.  A  quantitative  depiction  of  this  is 
seen  in  Figure  21. 
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Figure  21 .  Diffraction  of  the  optical  mode  compared  with  undulator  gap  width. 


The  straight  lines  in  the  top  half  of  Figure  21  depict  eleven  different  values  of  the 
beam  pipe  radius  (assumed  to  be  approximately  the  same  as  the  undulator  gap).  For  each 
undulator  gap,  the  color  scale  on  the  right  indicates  the  corresponding  wavelength.  The 
corresponding  optical  mode  for  each  wavelength  is  depicted  by  the  curved  lines  in  the 
lower  half  of  Figure  21,  showing  the  1/e  radius  of  the  Gaussian  optical  field  that  is  stored 
inside  a  concentric  optical  cavity.  The  maximum  distance  along  the  horizontal  axis  -390 
mm,  which  is  the  approximate  length  of  the  beam  pipe  that  will  be  inside  of  the 
undulator.  In  each  case,  the  1.8  mm  optical  mode  waist  radius  in  the  middle  of  the 
undulator  was  held  constant,  meaning  that  the  radius  of  curvature  of  the  mirrors  would 
change  for  each  case.  Power  loss  was  determined  by  calculating  the  power  clipped  from 
the  Gaussian  mode  by  the  inside  of  the  beam  pipe  at  the  ends  of  the  undulator.  As  can  be 
seen,  the  wide  28  mm  gap  (red  line)  gives  a  short  wavelength  (66  pm)  and  aids  in 
keeping  the  beam  pipe  far  from  the  optical  mode  to  cause  negligible  power  loss  due  to 
clipping.  Likewise  the  narrow  gap  (purple,  13  mm)  results  in  a  longer  wavelength  (180 
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^m),  which  leads  to  a  nearly  24  pereent  round  trip  loss  in  the  eavity  solely  due  to  elipping 
on  the  undulator  beam  pipe.  This  does  not  take  into  aeeount  other  possible  areas  of 
elipping  in  the  eavity  that  may  also  degrade  the  performanee.  To  maintain  a  round  trip 
loss  of  less  than  one  pereent  on  this  portion  of  the  beam  pipe,  the  longest  wavelength  that 
ean  be  aehieved  is  ~  126  pm  eorresponding  to  a  16  mm  gap.  As  the  gap  is  widened,  the 
losses  deerease,  but  so  does  the  gain,  sinee  the  deerease  in  transverse  motion  of  the 
eleetrons  also  lowers  the  eoupling  between  the  eleetrons  and  the  light.  This  leads  to  a 
total  effeetive  range  of  ~60  pm^l26  pm  when  just  eonsidering  elipping  from  the 
undulator  beam  pipe  [35]. 

F.  ELECTRON  BEAM  RADIUS 

To  maximize  energy  transfer,  the  eleetron  beam  and  optieal  mode  transverse 
profiles  must  properly  overlap  during  the  full  transit  of  the  undulator.  For  optimal 
bunehing,  gain,  and  extraetion,  the  eleetron  beam  should  be  well  inside  of  the  optieal 
mode.  Several  examples  using  GINGER  are  eonsidered  in  Figure  22. 
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Figure  22.  GINGER  simulations  of  electron  beam  (blue)  and  optical  mode  (red) 
transverse  profiles  for  2  mm,  1  mm,  and  0.55  mm  initial  electron  beam 
radius.  The  horizontal  axis  is  the  distance  z  along  the  undulator,  the 
vertical  axis  is  the  transverse  rms  dimension. 

Three  different  initial  electron  beam  radii  were  simulated  to  show  the  effect  of 
electron  beam  focusing  on  FEL  performance.  While  focusing  the  beam  aids  in 
maintaining  proper  overlap  inside  the  optical  mode,  over  focusing  can  be  detrimental.  A 
wider  initial  electron  beam  results  in  many  electrons  interacting  with  a  portion  of  the 
optical  field  with  lower  intensities,  reducing  bunching.  Also,  a  wider  initial  electron  beam 
requires  more  focusing  which  causes  a  large  spread  in  thereby  inducing  a  larger  spread 
in  initial  phase  velocities,  thus  inhibiting  the  bunching  mechanism.  In  the  top  graph  in 
Figure  22,  the  electron  beam’s  initial  radius  (blue)  is  2  mm  and  produces  a  weak  field 
gain  of  G~ll%.  In  the  middle  graph  the  electron  beam’s  initial  radius  is  1  mm,  which 
leads  to  an  increased  G=20%.  The  bottom  graph  is  the  nominal  case  that  is  being 
discussed  in  this  thesis  with  an  initial  electron  beam  radius  of  0.55  mm  at  the  beginning 
of  the  undulator,  which  leads  to  G~22%. 
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As  the  electron  beam  decreases  in  radius,  it  overlaps  the  portion  of  the  optical 
mode  with  higher  peak  fields.  This  increases  the  PEL  interaction  and  the  bunching 
mechanism,  thereby  improving  the  performance.  Figure  23  shows  simulation  results  for 
the  gain  for  an  initial  2  mm  radius  (first  case  in  Figure  22)  using  GINGER.  The 
simulation  was  started  with  an  initial  field  that  roughly  corresponds  to  spontaneous 
emission.  For  the  first  few  hundred  passes  the  simulation  shows  a  noisy  result  for  the  gain 
but  settles  down  to  ~10%  gain  for  the  next  few  hundred  passes.  After  -500  passes,  the 
gain  decreases  due  to  saturation  at  high  power.  When  the  gain  is  equal  to  the  round  trip 
losses  in  the  cavity  the  FEE  has  reached  steady-state  saturation. 


Figure  23.  GINGER  gain  as  a  function  of  pass  number  for  2  mm  initial  electron 

beam  radius. 

Figure  24  shows  simulation  results  for  the  gain  for  an  initial  1  mm  radius 
corresponding  to  the  second  case  in  Figure  22. 
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Figure  24.  GINGER  gain  as  a  function  of  pass  number  for  1  mm  initial  electron 

beam  radius. 
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For  the  1  mm  initial  radius  case,  saturation  is  reached  at  -'200  passes  due  to  the  larger 
-20%  weak  field  gain  caused  by  the  better  overlap  between  the  electron  beam  and  the 
optical  mode.  The  0.55  mm  initial  radius  case  has  even  higher  gain  and  therefore  reaches 
saturation  after  only  100  passes.  This  is  the  expected  result  in  the  experiment. 

G.  OPTICAL  MODE  SHAPE 

The  hole  out-coupling  technique  that  is  designed  for  the  optical  cavity  of  the 
Niowave  FEL  can  have  a  negative  effect  of  distorting  the  transverse  shape  of  the  optical 
mode.  Transferring  energy  away  from  the  fundamental  mode  can  decrease  the  coupling 
between  the  optical  mode  and  the  electron  beam,  reducing  the  FEL  gain.  The  size  of  the 
hole  needs  to  be  small  enough  so  the  round-trip  losses  do  not  exceed  the  weak  field  gain 
and  the  wavefronts  are  not  distorted  so  much  that  the  EEL  efficiency  decreases 
significantly.  However,  smaller  holes  produce  larger  diffraction  when  the  light  exits  the 
cavity,  which  must  also  be  considered.  Eigure  25  shows  the  transverse  electric  field 
profiles  at  the  undulator  entrance  for  a  range  of  hole  out-coupling  radii  from  0  mm  to  4 
mm  in  the  center  of  the  outcoupling  mirror. 
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Figure  25.  GINGER  simulation  results  for  transverse  eleetric  field  at  undulator 

exit  for  inereasing  hole  radii. 


Figure  25  depicts  the  transverse  electric  field  at  the  undulator  entrance  for  a  range 
of  output  coupler  hole  sizes.  The  side  lobes  at  larger  radii  (lower  right  in  Figure  25)  are 
due  to  the  energy  that  is  transferred  to  the  higher  order  modes.  Losses  and  gain  are 
calculated  in  the  weak  field  regime.  The  wavelength  that  was  simulated  was  the  nominal 
-100  pm  for  this  FEL  leading  to  a  radius  of  -27  mm  (1/e  in  the  field)  on  the  mirror.  The 
amount  of  gain  degradation  due  to  the  enlarged  hole  size  is  nearly  negligible  as  the  gain 
only  drops  -1%  at  the  largest  radius,  but  the  cavity  losses  are  large  (-15%).  Reiterating, 
GINGER  simulations  do  not  include  pulse  effect  so  the  actual  gain  will  be  lower  and 
could  be  below  this  loss  for  the  large  hole  radii.  For  this  design  a  hold  radius  of  0.5  mm 
will  be  initially  used  because  it  maintains  low  losses.  After  initial  operation,  the  hole 
radius  can  be  enlarged  to  increase  THz  output. 
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H.  AVERAGE  OUTPUT  POWER 


At  saturation,  the  amount  of  power  that  is  lost  by  the  eleetron  beam  is  equal  to  the 
output  power  of  the  PEL.  The  extraction,  which  is  the  fraction  of  the  electron  beam 
power  converted  to  light,  can  be  estimated  by  1/2V~5%  for  V=10  [30].  Because  of  the 
slippage  of  the  short  optical  and  electron  pulses,  this  should  be  considered  an 
overestimate.  All  the  codes  used  have  shown  an  extraction  of  '-'3-4%,  which  corresponds 
to  an  --2  kW  average  output  power. 

1.  OPTICAL  CAVITY  ALIGNMENT  TOLERANCES 

Ideally  the  optical  cavity  is  aligned  sufficiently  collinear  with  both  the  axis  of  the 
undulator  and  the  trajectory  of  the  electron  beam.  In  this  section,  the  tolerance  on  the 
optical  cavity  alignment  will  be  determined  by  simulating  the  FEE  interaction  with  either 
a  transverse  shift  or  angular  tilt  of  one  mirror.  Either  of  these  will  cause  a  misalignment 
of  the  electron  beam  and  the  optical  mode,  leading  to  decreased  interaction  length 
between  the  electrons  and  the  optical  mode,  thus  lowering  the  gain.  GINGER  was  not 
utilized  for  these  simulations  because  it  presumes  an  axisymmetric  radiation  field. 
Rather,  simulations  were  done  using  4DPULSEN  to  discern  the  effect  of  both  the 
transverse  shift  and  angular  tilt  offsets  (see  Figures  26  and  27). 
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Gain  dependence  on  mirror  shift 


Figure  26.  4DPULSEN  simulation  results  for  weak  field  gain  as  a  function  of 

shift  of  one  of  the  cavity  mirrors. 

Figure  26  shows  the  gain  plotted  against  the  transverse  mirror  offset,  where  the 
green  line  indicates  the  minimum  gain  threshold  of  six  percent  round  trip  loss  in  the 
cavity.  As  the  mirror  shift  increases,  the  gain  degrades  until,  at  -0.25  mm  offset,  the  gain 
is  below  threshold  and  the  FEL  will  not  work. 

This  same  analysis  is  done  for  the  angular  tilt  and  can  be  seen  in  Figure  27. 
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Gain  dependence  on  mirror  tilt 


Figure  27.  4DPULSEN  simulations  of  weak  field  gain  as  a  funetion  of  tilt  of 

one  of  the  eavity  mirrors. 

As  the  mirror  is  misaligned  by  an  angle  of  0  0.15  mrad,  the  gain  deereases  from  its 

maximum  value  of  ~13%  to  approximately  threshold.  This  gives  a  value  of  initial 
alignment  that  must  be  met  by  the  alignment  technique. 

As  before,  for  each  of  the  data  points  in  the  Figure  27  and  26,  the  amount  of 
desynchronism  has  been  optimized  to  maximize  the  gain.  The  change  in  desynchronism 
as  a  function  of  misalignment  can  be  seen  in  Figure  28. 


52 


Figure  28.  The  desynchronism  curves  for  different  mirror  shifts  (on  left)  and 

mirror  tilts  (on  right). 


As  the  cavity  mirror  is  misaligned,  the  maximum  gain  decreases  and  occurs  at  a 
smaller  desynchronism.  For  no  mirror  shift  (left  in  Figure  28)  Aym=0.0  mm  (blue),  the 
peak  gain  occurs  at  d=0.02  whereas  for  Aym=0.23  mm  (black)  the  peak  is  at  d=0.0\.  For 
no  mirror  tilt  (right  in  Figure  28),  A0m  -0.0  mrad  (blue),  the  maximum  gain  is  at  d=0.02, 
whereas  for  the  A0m=0.13  mrad  the  peak  occurs  at  £/=0.01.  This  change  in  d  (from  0.01  to 
0.02)  corresponds  to  a  4.5  pm  shift  in  the  cavity  length.  Knowing  the  mirror  shift  and  tilt 
tolerances  enables  the  engineering  design  of  an  alignment  procedure,  which  will  be 
described  in  a  later  chapter. 
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IV.  ELECTRON  BEAM  PHYSICS 


A.  ELECTRON  BEAM  QUALITY 

The  electron  beam  is  made  up  microbunches  of  free  electrons  that  are  co- 
propagating  with  each  other.  Each  electron  has  a  position  r  and  a  velocity  r  with  respect 
to  the  beam  centroid,  which  describes  its  position  and  momentum.  In  the  lab  reference 
frame,  the  values  of  f  are  seen  as  a  distribution  of  angles  with  respect  to  the  propagation 
vector  of  the  beam  centroid. 

1.  Transverse  Emittance 

Emittance  is  defined  as  the  area  of  the  electron  beam  in  phase  space  and  is  given 
by  f  =  (Tr(Jg,  where  a,-  is  the  rms  transverse  spread  in  positions  and  ae  is  the  rms  spread  in 
angles  from  the  propagation  direction.  Rather  than  looking  at  the  particles  in  position 
space  (where  there  3  dimensional  location  is  considered),  it  is  often  better  to  analyze  the 
beam  in  phase  space,  or  where  a  particles  momentum  and  position  are  compared.  A 
typical  phase  space  diagram  can  be  seen  in  Eigure  29. 


Eigure  29.  Typical  phase  space  diagram  of  electron  bunch  where  x  is  the 

particles  transverse  position  x  ’  is  the  particles  momentum  in  that 
direction,  from  [36]. 
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In  Figure  29  the  eleetron  buneh  has  a  spread  in  both  angles  and  positions  (in  the 
x-x’  plane).  By  Liousville’s  theorem,  the  volume  in  phase  space  is  constant  in  the  bunch 
frame.  Emittance  is  often  measured  at  low  energies,  but  at  the  undulator  entrance  (where 
the  necessary  beam  quality  is  usually  defined),  the  electron  beam  is  at  a  much  higher 
energy.  Due  to  relativistic  effects,  the  emittance  is  dependent  on  the  speed  of  the 
electrons.  The  normalized  emittance,  removes  the  dependence  on  the 

velocity  and  supplies  a  measure  of  beam  quality  that  is  constant  throughout  a  whole  beam 
line  (as  long  as  the  number  of  particles  lost  is  inconsequential). 

The  phase  space  volume  remains  constant  as  an  electron  bunch  transits  a  beam 
line  but  may  change  shape.  A  typical  beam  line  component  is  a  magnet  to  supply 
focusing  to  the  electron  beam.  Figure  30  depicts  the  evolution  of  the  phase  space  as  a 
beam  is  focused  (in  the  x  transverse  plane)  through  a  waist. 

Beam  envelope  along  beamline 


Figure  30.  Electron  beam  phase  space  evolution  as  it  passes  through  a  beam 

waist,  from  [36]. 


In  the  left  side  of  Figure  30,  the  beam  is  converging  towards  the  waist.  This  can  be  seen 
in  the  phase  space  diagram  in  the  lower  left,  where  electrons  with  negative  positions 

(x<0)  have  positive  angles  (x’>0)  and  the  electrons  with  positive  positions  have  negative 

56 


angles.  As  the  beam  propagates  to  the  waist,  angular  spread  remains  constant  while  the 
maximum  radius  decreases;  this  induces  an  effective  sheer  on  the  phase  space  area.  At 
the  waist,  the  radius  is  at  a  minimum,  as  can  be  seen  by  the  reduced  transverse  extent  of 
the  phase  space  diagram  in  the  center  of  the  figure.  But,  the  beam  has  both  positive  and 
negative  angles.  As  the  beam  expands  beyond  the  waist,  the  electrons  with  positive 
positions  (x>0)  now  have  positive  angles  (x’>0),  and  those  with  negative  positions  have 
negative  angles,  increasing  the  transverse  size  of  the  electron  beam.  By  Liousville’s 
theorem  the  area  remains  constant.  A  large  emittance,  or  large  area  in  phase  space, 
inhibits  the  ability  to  focus  the  beam  to  a  small  radius. 

For  the  FEL  to  work  optimally,  the  electron  beam  should  be  well  inside  the 
optical  mode  along  the  length  of  the  undulator.  This  places  geometrical  constraints  on  the 
electron  beam  waist  radius  and  angular  spread,  and  therefore  the  emittance.  Applying 
these  concepts  leads  to  a  characteristic  emittance  £c  ~  yXI\n  above  which  weak  field 
gain  is  diminished  [37].  For  the  Niowave  FEL  parameters,  £c~135  mm-mrad  and  is 
several  times  larger  than  the  nominal  28  mm-mrad  used  in  the  simulations.  The 
simulation  code  4DPULSEN  confirms  that  an  emittance  £„=100  mm-mrad  decreases 
weak  field  gain  from  approximately  14  percent  down  to  approximately  10  percent  and  the 
EEL  fails  to  operate  when  emittance  is  increased  to  approximately  200  mm-mrad. 

2.  Energy  Spread 

During  the  generation  and  acceleration  of  the  electron  beam,  different  electrons 
will  gain  different  amounts  of  kinetic  energy  due  to  entering  the  accelerating  structures  at 
slightly  different  RE  phases.  In  addition  to  affecting  EEL  performance  (as  was  previously 
described  in  Chapter  III),  this  energy  spread  can  cause  problems  with  beam  transport.  As 
magnetic  fields  along  the  beam  line  are  used  to  bend  and  focus  the  electrons,  a  spread  in 
energy  leads  to  individual  electrons  being  affected  differently,  since  the  magnetic  force 
depends  on  their  energy.  These  effects  and  their  mitigation  will  be  discussed  in  later 
sections. 
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3. 


Desired  Electron  Beam  Quality 


The  electron  beam  quality  requirements  for  the  Niowave  THz  PEL  are  relatively 
lax  due  to  the  long  wavelength  of  the  FEE  and  the  short  undulator.  The  longer 
wavelength  allows  for  a  larger  energy  spread  as  bunching  on  the  order  of  a  wavelength  of 
light  is  easier.  A  shorter  undulator  increases  the  gain  spectrum  bandwidth,  -'1/A,  so  the 
FEE  has  less  sensitivity  to  energy  spread  from  the  accelerator.  Also,  with  such  a  short 
undulator,  a  beam  with  a  large  emittance  can  still  be  focused  properly  through  it  without 
scraping  on  the  beam  pipe.  The  electron  beam  delivered  to  the  undulator  must  be 
matched  to  fit  inside  the  optical  mode  to  supply  proper  overlap  along  the  undulator,  as 
was  described  in  the  Chapter  III.  Simulation  results  indicate  that  the  Niowave  FEE  can 
tolerate  a  large  transverse  emittance  of -100  mm-mrad  with  up  to  -1.7%  energy  spread. 
For  the  nominal  emittance  value  of  28  mm-mrad,  an  energy  spread  up  to  three  percent  is 
tolerable. 

B.  FOCUSING:  QUADRUPOLE  MAGNETS 

In  order  to  focus  the  electron  beam,  quadrupole  magnet  is  formed  by  four 
hyperbolic  pole  face  electromagnets  arranged  symmetrically  around  the  beam  pipe  with 
alternating  polarity  magnets.  A  depiction  of  the  arrangement  can  be  seen  in  Figure  31. 
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Figure  3 1 .  Depiction  of  magnetic  field  of  quadrupole  magnet,  after  [38]. 

The  magnets  are  placed  axially  around  the  beam  line  creating  a  magnetic  field  at  the  site 
of  the  quadrupole  of  the  form 

Eq.  4.1 

Uq 

By  =  Bq-,  Eq.  4.2 

dq 

where  Bg  is  the  peak  magnetic  field  next  to  the  pole  face,  and  (x,y)  are  measured  from  the 
axis,  and  Ug  is  the  distance  from  the  center  axis  to  the  magnet  face.  The  strength  of  the 
focusing  is  controlled  by  the  combination  of  the  interaction  length  of  the  quadrupole  as 
well  as  the  peak  magnetic  field.  In  Figure  31,  the  electron  beam  is  going  into  the  page,  so 
this  magnet  would  supply  focusing  in  the  y  direction  and  defocusing  in  the  x  direction. 
Due  to  the  focusing  /  defocusing  nature  of  quadrupole  magnets,  they  are  commonly 
placed  in  groupings  of  two  or  three  magnets  to  supply  the  necessary  focusing  in  both 
planes.  One  example  of  a  quadrupole  “triplet”  can  be  seen  in  Figure  32. 
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Figure  32.  Quadrupole  triplet  for  transverse  beam  focusing,  from  [23]. 


The  strength  of  the  magnets  and  the  size  of  the  drift  space  between  quadrupoles 
must  be  optimized  to  supply  the  necessary  focusing  for  the  beam  at  the  designed  focal 
point  later  in  the  beam  line.  With  each  quadrupole  acting  as  a  lens  with  focal  lengths  fyji, 
^3,  respectively,  it  can  easily  be  shown  that  the  focal  length  of  three  quadrupoles  in  series 
is 


_ /1/2/3 _ 

/3(/2+/i~'5i)+/i/2+S2/2+S2/i-SiS2  ’ 


Eq.  4.3 


where  and  S2  are  the  drift  spaces  between  subsequent  quadrupoles  [39].  These 
quadrupoles  will  be  used  to  help  focus  the  electron  beam  into  the  undulator  of  the  FEE  to 
create  the  correct  overlap  between  the  optical  mode  and  the  electron  beam  and  create  the 
minimum  radius  needed. 


C.  BENDING:  DIPOLE  MAGNET  DESIGN 

A  dipole  magnet  supplies  a  constant  magnetic  field  across  the  beam  path  to  bend 
the  electron  beam. 
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1.  Dipole  Design 


One  design  of  a  dipole  magnet  eonsists  of  an  iron  yoke  wrapped  with  coils  of 
wire;  when  current  is  applied  through  the  coils  it  produces  an  approximate  uniform 
magnetic  field  across  a  gap.  The  size  and  shape  of  the  iron  yoke  controls  the  length, 
width,  and  gap  size  of  the  dipole.  The  current  through  the  wire  controls  the  magnetic  field 
strength  and  polarity.  The  magnetic  field  is  enhanced  by  the  permeability  of  the  iron  yoke 
[40]. 

A  cross  section  of  two  typical  electromagnetic  dipoles  can  be  seen  in  Figure  33. 


Figure  33.  Cross-sectional  view  of  dipole  magnets,  from  [41]. 


On  the  left  is  a  “C-core”  magnet  that  has  one  open  side  of  iron;  on  the  right  is  a  “H-type” 
yoke,  which  is  symmetrical  about  the  center  axis.  The  coils  of  wire  are  denoted  by  the 
boxes  surrounding  the  solid  yokes  and  the  magnetic  field  directions  and  magnitudes  are 
given  by  the  direction  and  size  of  the  arrows.  Inside  the  yokes,  the  color  denotes  the 
magnetic  field  strength. 

2.  Simple  Dipole  Physics 

The  geometry  of  the  interaction  between  an  electron  and  a  uniform  magnetic  field 
can  be  seen  in  Figure  34. 
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n  n  n  n  n 


Figure  34.  Depiction  of  electrons  motion  through  a  uniform  magnetic  field, 

after  [42], 


The  electron  enters  from  the  left  where  it  begins  to  interact  with  a  uniform  magnetic  field 
of  magnitude  Bd  (into  the  page)  for  a  length  Ld.  The  relativistic  Lorentz  force  alter  the 
electron’s  trajectory  by  an  angle  0.  The  relativistic  Lorentz  force  is 

—  (ymv)  =  — ev  x  ,  Eq.  4.4 


leading  to 


dVx  ^  ,  dVz  ^  eBd  ^  ^ 

-  X  H - Z  =  -  (WyX  —  VyZ)  . 

dt  dt  ym  ^  ^  ^ 


Eq.  4.5 


Differentiating  again  with  respect  to  time  gives 


d^Vx 

dt^ 


I 

X  + 


rf  Vz 
dt^ 


Eq.  4.6 


which  can  be  separated  into  two  differential  equations,  each  in  the  form  of  a  simple 
harmonic  oscillator  with  frequency  0)^  =  eB^i/yvn.  The  solution  can  thus  be  written 

V;,  = -)gcsin(mdt),  Eq.  4.7 

V2  =  cos(a)at),  Eq.  4.8 

where  the  initial  conditions  are  Vzo=Ac,  Vxo=0.  The  deflection  angle  is  given  by  tan(0)  = 
Vx/Vy  =  tan(a)dt),  leading  to 

=  Eq.  4.9 
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The  deflection  angle  depends  on  the  strength  of  the  magnetic  field  {Bd),  the  electron’s 
kinetic  energy  (y)  and  the  time  (t)  spent  in  the  field  (i.e.,  the  length  of  the  dipole).  While 
the  energy  of  the  electron  is  set  by  the  design  of  the  particle  accelerator,  the  length  of  the 
dipole  and  the  strength  of  the  magnetic  field  must  be  controlled  to  induce  the  proper 
change  in  trajectory. 

3.  Fringe  Fields 

The  entire  derivation  leading  to  Eq.  4.9  is  based  on  an  idealistic  uniform 
rectangular  magnetic  field  with  “hard  edges.”  Realistically,  dipoles  have  fringe  fields  that 
extend  past  the  edges  of  the  magnets  into  the  surrounding  area.  These  can  be  seen  in 
Figure  33  on  both  types  of  magnets.  This  effectively  increases  the  dipole  length  resulting 
in  an  increased  deflection  angle.  This  effect  must  be  taken  into  account  when  designing 
the  strength  of  the  magnet. 

4.  Energy  Spread 

A  realistic  beam  is  comprised  of  billions  of  electrons  with  spreads  in  position, 
angle,  and  energy.  An  electron’s  position  and  angle  with  respect  to  the  beam  centroid  will 
be  translated  through  the  dipole  without  affecting  the  beam’s  properties,  but  the  energy 
spread  in  the  electron  beam  will  cause  an  enlarged  spread  in  angles  after  the  dipole  since 
the  energy  of  an  electron  affects  its  deflection  angle  according  to  Eq.  4.9.  An  example  of 
this  spread  can  be  seen  in  Figure  35. 
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Dipole  magnet 


Figure  35.  Electron  beam  dispersion  from  a  dipole  magnet  due  to  energy 

spread. 

Figure  35  depicts  a  2  mm  radius,  8  MeV  electron  beam  with  a  -1.7%  percent  energy 
spread  entering  the  dipole  (green)  from  the  left.  After  the  dipole  the  electron  trajectories 
are  altered,  but  due  to  the  enlarged  spread  in  angles  the  transverse  size  of  the  beam 
grows.  This  can  be  a  problem  when  trying  to  focus  the  electron  beam  inside  the 
undulator. 

D.  COMBINED  FUNCTION  MAGNETS 

Some  magnets  can  combine  multiple  functions  (i.e.,  focusing  and  bending)  by 
virtue  of  their  shape.  In  the  case  of  the  Niowave  TFlz  FEE,  a  single  magnet  can  be  used  to 
bend  the  electron  beam  while  controlling  the  dispersion  caused  by  the  energy  spread. 
Two  techniques  can  used  to  control  the  angular  deflection  of  electrons  as  they  pass 
through  a  magnet:  either  the  time  the  electron  interacts  with  the  magnetic  field  can  be 
preferentially  changed,  or  the  strength  of  the  magnetic  field  can  be  altered  in  such  a 
manner  that  different  energy  electrons  will  interact  with  different  magnitude  magnetic 
fields. 
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1. 


Tilted  End  Faces 


To  alter  the  interaction  length  the  individual  electrons  have  with  the  magnetic 
field,  the  shape  of  the  magnetic  field  region  needs  to  be  adjusted  from  the  simple 
rectangular  area  that  was  previously  considered  (see  Figure  36). 


Figure  36.  Dipole  magnet  with  rotated  pole  faces. 

For  example,  Figure  36  depicts  a  top  view  of  a  dipole  magnet  designed  to  alter  the 
trajectory  of  an  electron  beam  by  an  angle  a.  In  this  geometry  the  nominal  end  face 
directions  (dotted  line)  are  orthogonal  to  the  propagation  direction  of  the  electron  beam. 
In  both  the  left  and  right  diagram  the  entrance  and  exit  faces  of  the  magnet  are  tilted  from 
this  nominal  end  face  by  an  amount  cp^  and  (pout,  respectively.  When  an  electron  beam 
enters  the  dipole,  the  higher  energy  electrons  will  be  less  affected  by  the  magnetic  field 
and  will  tend  to  go  towards  the  outside  of  the  curve.  To  compensate  for  this,  the  end  faces 
are  tilted  in  such  a  manner  so  that  the  higher  energy  electrons  will  interact  with  the  dipole 
longer  than  the  lower  energy  electrons.  This  supplies  a  transverse  focusing  ability  to  the 
dipole  and  will  both  bend  the  electron  beam  as  well  as  focusing  it  in  the  plane  of  the  turn, 
counteracting  the  dispersion. 

2.  Canted  Pole  Faces 

Another  method  to  preferentially  control  the  angular  deflection  of  the  individual 


electrons  is  to  cant  the  transverse  pole  faces.  This  involves  turning  the  pole  faces  of  the 
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magnets,  which  are  normally  parallel  to  each  other,  thus  adding  a  transverse  dependence 
to  the  strength  of  the  magnetic  field,  as  seen  in  Figure  37. 


Figure  37.  End  view  of  canted  dipole,  from  [43]. 

The  end  view  of  a  canted  dipole  in  Figure  37  shows  the  angle  of  the  pole  faces  with 
respect  to  the  direction  of  travel.  Much  like  the  tilted  end  face  design,  this  magnet 
preferentially  induces  more,  or  less,  force  on  an  electron  depending  on  its  energy.  Rather 
than  altering  the  amount  of  time  the  electron  travels  through  the  dipole  based  on  its 
energy,  this  magnet  alters  the  strength  of  the  magnetic  field  to  adjust  the  output  angle 
accordingly. 

E.  NIOWAVE  BEAMLINE  GEOMETRY 

The  electron  beam  needs  to  be  transported  from  the  final  accelerating  structure  of 
the  linac  to  the  undulator  and  must  be  focused  in  the  correct  manner  necessary  for  the 
FEE  to  work.  The  geometry  of  the  Niowave  beamline  after  the  linac  can  be  seen  in 
Eigure  38. 
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Side  View  3m  0.6m  0.81m  0.71m 


Figure  38.  Side  and  top  profde  view  of  Niowave  beamline. 

The  electron  beam  enters  from  the  left  and  transits  three  meters  of  beam  pipes  before 
turning  through  two  horizontal  23°  bends  (one  +23°  and  the  other  -23°)  to  be  injected 
into  the  undulator.  After  transiting  through  the  undulator,  it  then  goes  through  a  vertical 
45°  bend  downward  towards  the  ground,  where  the  electrons  are  sent  into  a  beam  dump. 
In  the  middle  of  the  undulator,  the  electron  beam  must  be  focused  to  a  nominal  0.45mm 
waist  radius  (yrms~Xrms~0.32  mm).  To  accomplish  this,  both  quadrupole  and  dipole 
magnets  will  be  utilized. 

F.  BEAMLINE  MODELING:  GENERAL  PARTICLE  TRACER 

A  beamline  model  can  be  used  to  design  the  necessary  components  to  deliver  an 
electron  beam  to  the  undulator  with  the  proper  qualities. 

I.  Code  Description 

A  commercial  particle  tracking  code  is  used  to  design  the  beamline  and  quantify 
the  parameters  of  the  beamline  components.  General  Particle  Tracer  is  a  complete 
software  package  produced  by  Pulsar  Physics  to  study  charged  particle  dynamics  in 
electromagnetic  fields  [44].  An  input  file  describes  the  beamline  components’  positions 
and  strengths,  the  electron  beam’s  initial  properties  (energy  spread,  position  spread, 
emittance,  etc...),  as  well  as  the  diagnostic  locations  along  the  beamline  (and  what 
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information  is  to  be  outputted).  It  maps  out  the  entire  three-dimensional  eleetromagnetie 
field  and  then,  using  a  fifth  order  Runge-Kutta  method,  calculates  the  individual  particle 
trajectories  through  the  beam  line. 

2.  Optimization  with  Only  Quadrupoles 

To  supply  the  necessary  focusing  to  the  electron  beam,  a  set  of  three  quadrupole 
magnets  are  utilized.  As  seen  in  Figure  38,  the  distance  from  the  start  of  the  simulation  at 
the  linac  exit  (z=0)  to  the  middle  of  the  undulator  is  ~4.4  m.  The  quadrupole  triplet  will 
be  placed  before  the  first  dipole.  GPT  contains  built  in  functions  to  optimize  parameters 
to  focus  both  transverse  directions  to  the  necessary  spot  size  and  location  by  comparing 
various  quadrupole  strengths  and  positions.  The  electron  beam  parameters  that  were 
simulated  are  listed  in  Table  5. 


Beam  kinetic  energy 

8MeV 

Initial  electron  beam  radius  (Xrms,Trms) 

2  mm 

Normalized  emittance,  rms 

28  mm-mrad 

Energy  spread  (Ay/y) 

1.66% 

Tables.  GPT  simulation  parameters. 


This  resultant  profile  from  optimizing  the  quadrupole  parameters  to  focus  in  the 
middle  of  the  undulator  can  be  seen  in  Figure  39. 
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Figure  39.  Electron  beam  radius  (rms)  through  beamline. 

A  reasonable  assumption  for  the  electron  beam’s  initial  radius  after  the  linac  is  taken  to 
be  2  mm.  For  this  simulation  GPT  was  tasked  with  minimizing  the  beam  radius  in  both 
transverse  directions  at  z=4.4  m  by  adjusting  quadrupole  parameters.  The  quadrupole 
parameters  that  produced  this  result  are  listed  in  Table  6. 


F'  quadrupole  location 

0.5  m 

F‘  quadrupole  peak  magnetic  field  gradient 

1.72  T/m 

2"‘‘  quadrupole  location 

1.8m 

2"“  quadrupole  peak  magnetic  field  gradient 

-1.57  T/m 

quadrupole  location 

2.25  m 

3'^‘‘  quadrupole  peak  magnetic  field  gradient 

1.64  T/m 

Table  6.  Optimized  quadrupole  parameters  with  the  distances  measured 

from  the  end  of  the  linac. 
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These  parameters  lead  to  a  beam  size  of  0.88  mm  and  0.57  mm  in  the  x  and  y  direetion 
respeetively  at  4.4  m.  Both  are  larger  than  the  -0.32  mm  mentioned  previously  but  are 
smaller  than  the  1.3  mm  optical  mode  (components).  As  long  as  the  electron  beam  is  well 
inside  the  optical  mode,  the  PEL  interaction  works  as  desired.  The  quadrupoles  were 
assumed  to  have  a  length  of  3  cm,  which  is  a  reasonable  assumption  [45].  GPT  does  not 
include  fringe  field  effects  in  the  quadrupole  magnets  so  the  actual  peak  magnetic  field 
would  be  slightly  less  than  listed  in  Table  6. 

3.  Combining  Quadrupoles  and  Rectangular  Dipole  Magnets 

Adding  the  two  dipoles,  which  are  needed  to  bend  the  electron  beam  around  the 
cavity  mirrors,  will  cause  a  decrease  in  focusing  in  the  x-plane  due  to  dispersion.  The  two 
dipoles  are  separated  by  -54  cm  and  the  magnetic  fields  will  need  to  be  equal  strength 
but  opposite  polarity  to  supply  the  -0.23  m  offset  without  altering  the  final  propagation 
direction  into  the  undulator.  The  dipoles  were  assumed  to  be  20  cm  long  which  is  a 
reasonable  design  given  the  low  energy  of  the  electrons.  Without  the  addition  of  fringe 
field  effects  the  necessary  magnetic  field  strength  can  be  calculated  from  Eq.  4.7  and  is 
equal  to  -0.056  Tesla.  Adding  the  two  dipoles  leads  to  the  beam  path  seen  in  Eigure  40 
from  GPT. 


70 


Figure  40.  Electron  beam  trajectory  in  x  (top)  and  y  (bottom)  plane. 


This  dipole  chicane  supplies  the  required  translational  offset  but  the  dispersion 
leads  to  a  beam  size  of  ~4  mm  by  -“0.55  mm  in  the  x  and  y  directions  respectively.  Note 
that  the  transverse  size  in  y  is  not  altered  by  the  dipoles,  as  expected.  This  would  cause 
approximately  half  of  the  electron  beam  to  be  outside  the  optical  mode  inside  the 
undulator,  which  would  severely  degrade  FEE  performance.  The  effect  the  dipole  has  can 
be  seen  in  a  graph  of  the  electron  beam’s  radius  in  the  x  direction  in  Figure  41. 
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Quadrupole  Dipole  Middle  of 


Figure  41 .  Transverse  eleetron  beam  size  along  beam  path. 

The  electron  beam  size  follows  the  same  trend  as  seen  in  Figure  39  until  the  third  meter, 
where  the  radius  starts  to  increase  due  to  the  dipole  dispersion.  As  the  beam  transits  the 
second  dipole,  the  opposite  polarity  field  induces  another  spread  in  deflection  angles,  but 
this  one  is  equal  and  opposite  of  the  first  cancelling  out  the  spread  in  angles.  It  is  in  the 
drift  space  between  the  two  dipoles  that  the  electron  beam’s  transverse  size  is  increased. 
Due  to  this,  the  distance  was  minimized  for  the  given  angle  and  geometry. 

4.  Quadrupoles  and  Combination  Magnets 

As  was  described  in  section  D.2,  the  use  of  combination  magnets  can  be  used  to 
counteract  the  dispersion  created  by  the  dipole.  For  this  type  of  magnet,  the  nominal 
length  of  the  dipole  and  the  necessary  bending  angle  describes  the  radius  of  curvature. 
GPT  contains  a  built  in  model  of  this  type  of  magnet;  the  geometry  can  be  seen  in  Figure 
42. 
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Figure  42.  Geometry  of  eombination  magnet,  from  [44], 

In  Figure  42,  the  electron  beam  enters  from  the  left,  interacts  with  the  magnetic  field  and 
is  turned  through  a  circle  of  radius  R  before  exiting  to  the  right.  Nominally  the  pole  faces 
are  orthogonal  to  the  incident  electron  beam  direction  and  the  exiting  electron  beam 
direction.  To  supply  the  necessary  focusing  these  pole  faces  can  be  rotated  to  alter  the 
path  length  for  different  energy  electrons  undergo.  By  altering  this  angle  while  keeping 
the  dipole  placement  the  same  as  before  leads  to  a  properly  focused  beam  through  the 
undulator.  The  resulting  particle  trajectories  can  be  seen  in  Figure  43. 
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Figure  43.  Electron  beam  trajectories  in  x  (top)  and  y  (bottom)  plane  using 

combination  magnets. 


The  addition  of  the  combination  magnets  allows  the  electron  to  be  focused  in  the 
x-plane  into  the  undulator.  To  minimize  the  electron  beam  radius  in  the  middle  of  the 
undulator  in  this  case  required  that  both  end  faces  of  the  first  dipole  be  turned  15°  and  the 
end  faces  of  the  second  dipole  rotated  -10°  (as  described  in  Figure  42).  The  nominal 
length  of  the  dipole  is  still  20  cm  leading  to  the  same  magnetic  field  strength  of  -0.056 
Tesla  that  was  calculated  for  the  rectangular  dipole  magnets.  Fringe  fields  were  not 
considered  in  this  case  so  a  slightly  reduced  magnetic  field  would  be  needed  to  produce 
the  necessary  bending  angle.  With  this  combination  a  beam  size  of -1  mm  and  -0.55  mm 
in  both  X  and  y  directions  is  achieved  in  the  middle  of  the  undulator.  While  this  is  greater 
than  the  initial  design  goal,  it  is  well  inside  the  1.8  mm  optical  mode  radius  that  is 
designed  for  the  optical  cavity,  and  the  FEE  interaction  will  work  as  designed. 

After  the  FEE  interaction,  the  electron  beam  will  have  a  drastically  increased 
induced  energy  spread  due  to  the  relatively  high  extraction  in  the  FEE.  This  induced 
energy  spread  causes  the  beam  to  be  difficult  to  control  in  a  precise  manner  as  the 
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previous  chromatic  dispersion  problems  are  drastically  worsened  for  the  dipole  magnet 
bending  the  beam  towards  the  beam  dump.  The  results  from  a  GPT  simulation  of  the 
electron  beam  after  the  undulator  with  the  large  induced  energy  spread  can  be  seen  in 
Figure  44. 


Figure  44.  Trajectories  of  electrons  from  middle  of  undulator  through  last 

dipole. 

The  simulation  shows  the  electron  trajectories  in  the  vertical  plane  as  the  electrons 
would  be  removed  from  the  optical  cavity.  At  this  point  the  electron  beam  needs  to  be 
discarded  into  a  beam  dump.  Combination  magnets  can  be  used  to  attempt  to  mitigate  the 
dispersion  but  the  simplest  solution,  and  the  one  that  is  recommended  for  this  case,  is  to 
place  the  beam  dump  as  close  as  possible  after  the  last  dipole.  If  the  beam  dump  is  placed 
20  cm  after  the  last  dipole  the  electrons  beam  transverse  diameter  will  be  only  ~6  cm  by 
-“1.6  cm,  which  is  well  inside  the  four  inch  diameter  beam  pipe  that  is  being  used  already 
for  the  beam  pipe.  The  beam  spread  can  be  useful  to  prevent  damage  to  the  beam  dump 
from  the  focused  high  power  of  the  beam. 

a.  GINGER  Results 

The  electron  beam  parameters  from  the  GPT  simulations  can  be  used  as  input  into 
GINGER  to  estimate  the  performance  of  the  FEL  with  the  preceding  beam  line.  This 
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involved  calculating  the  electron  beam’s  size  and  angular  characteristics  at  the  beginning 
of  the  undulator  in  GPT  and  transferring  those  parameters  into  GINGER.  The  emittance 
of  a  beam  can  be  described  by  a  relation  with  the  emittance,  e  —  +  2ajxx'  + 

where  aj,  Pt,  Jt  are  the  Twiss  parameters  [36].  With  77^  =  (1  +  the 

twiss  parameters  at  the  undulator  entrance  created  by  the  beamline  components  from 
GPT  are  listed  in  Table  7. 


x(rms) 

y(rms) 

O-Tx 

ary 

Ptx 

hy 

0.81  mm 

0.87  mm 

1.27 

1.5 

0.26  m 

0.47  m 

Table  7.  GPT  derived  Twiss  parameters  of  electron  beam  at  undulator 

entrance. 


The  simulation  results  of  the  electron  beam  and  optical  mode  profiles  that  GINGER 
produces  can  be  seen  in  Figure  45. 


Figure  45.  Electron  beam  (red)  and  optical  mode  (blue)  transverse  profiles  for 
GINGER  simulations  using  the  electron  beam  parameters  from  GPT. 

The  non-symmetric  nature  of  the  electron  beam  (blue)  focusing  in  both  planes  is  shown 
but  the  beam  is  still  well  within  the  optical  field  (red).  GINGER  includes  betatron 
focusing  from  the  undulator,  which  decreases  the  electron  beam  waist  size  slightly  [46]. 
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This  provides  an  rms  waist  size  of  -0.24  mm  and  0.27  mm  in  the  x  and  y  direction 
respectively.  The  resulting  PEL  performance  can  be  seen  in  Figure  46. 


Gain 
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Figure  46.  Ginger  simulations  for  gain  for  the  Niowave  THz  FEE  using  electron 

beam  properties  from  GPT. 

These  electron  beam  focusing  parameters  closely  match  the  nominal  results  shown  in 
Figure  46  and  supply  a  weak  field  gain  of  -22%  in  approximately  the  same  number  of 
passes  as  the  nominal  case. 
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V.  HARDWARE  DESIGN 


A.  LINAC  COMPONENTS 

Niowave  Inc.  has  designed  and  built  linacs  for  various  U.S.  Navy  projects  in  the 
past.  One  project  has  incorporated  a  normal  conducting  copper  cavity  injector  using  a 
thermionic  cathode  and  a  superconducting  spoke  cavity  for  the  accelerating  structure, 
such  as  the  one  designed  for  the  Niowave  THz  PEL,  depicted  in  Figure  47. 


Figure  47.  CAD  drawing  of  linac  for  Niowave  THz  FEE. 

The  electron  beam  is  generated  in  the  copper  injector  (left)  before  entering  the 
superconducting  spoke  cavity  accelerator  (middle).  The  cutaway  shows  the  dewers  that 
surround  the  cavity  to  hold  the  outer  tank  of  liquid  nitrogen  and  the  inner  tank  of  helium 
to  maintain  the  proper  temperature  of  the  niobium.  The  electron  beam  exits  to  the  right 
after  being  accelerated  to  the  desired  8  MeV  energy.  It  is  after  this  point  that  the 
geometry  and  beamline  components  are  designed  in  previous  sections  to  transport  the 
electron  beam  in  a  controlled  way  to  the  undulator. 
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B.  TRANSPORT  BEAMLINE  COMPONENTS 


The  electron  beam  must  be  transported  from  the  linac  to  the  undulator  under 
vacuum,  bending  around  the  cavity  mirrors.  The  beam  must  also  be  focused  into  the 
undulator  to  properly  overlap  the  optical  pulse.  A  design  of  the  transport  beamline  for  the 
Niowave  THz  PEL  can  be  seen  in  Figure  48. 


Transport  Beamline 

- 0 - 

Top  View  ^  Quadruples 


Dipoles 


Undulator 


Figure  48.  Niowave  TPlz  FEE  transfer  beamline  highlighted  in  red. 


The  transport  beamline  consists  of  three  quadrupoles  that  supply  focusing  and  two 
combination  dipole  magnets  to  supply  bending  and  additional  focusing  into  the  undulator. 

1.  Beam  Pipes 

The  beam  pipes  that  transport  the  electron  beam  from  the  linac  to  the  optical 
cavity  must  be  large  enough  to  avoid  losses  of  the  beam  current  due  to  scraping,  while 
not  being  unreasonably  wide  as  to  increase  the  cost  of  the  beamline  magnetic 
components.  A  larger  beam  pipe  requires  a  stronger  magnet  due  to  the  larger  separation 
distances  involved.  As  can  be  seen  in  Figure  41,  the  largest  root-mean-square  (rms) 
radius  the  electron  bunch  achieves  in  either  transverse  direction  during  the  transport  to 
the  undulator  is  only  --7  mm.  To  allow  for  a  margin  of  error,  a  beam  pipe  diameter  of  --38 
mm  will  be  utilized  so  that  the  inside  pipe  wall  is  approximately  five  e-foldings  away 
from  the  beam.  It  is  estimated  that  only  ~0.1%  of  the  beam  will  be  lost  in  this  section  of 
the  beamline.  The  beam  pipe  will  need  to  be  custom  made  to  supply  the  necessary  -14° 
bend  while  also  connecting  with  the  optical  cavity  beam  pipe.  An  example  of  the  type  of 
custom  beam  pipe  that  will  be  necessary  is  shown  in  Figure  49. 
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Figure  49.  Diagram  of  necessary  custom  beam  pipe  with  14  degree  bends. 


The  two  bends  along  with  the  transition  from  the  transfer  beam  pipe  diameter  to  the 
optical  cavity  beam  pipe  diameter  may  need  to  be  custom  fabricated.  The  cost  of 
beampipes  is  on  the  order  of  a  few  thousand  dollars. 

2.  Quadrupoles 

The  quadrupoles  were  designed  to  supply  the  focusing  in  the  y  direction  while 
aiding  the  dipoles  in  focusing  in  the  x  direction.  The  three  quadrupoles  were  each 
assumed  to  be  3  cm  long  with  peak  magnetic  field  gradients  of  1.68  T/m,  -1.57  T/m,  and 
1.61  T/m  respectively.  For  permanent  magnet  quadrupoles,  a  pole  tip  field  of  1  T  is 
achievable,  which  in  conjunction  with  the  38mm  diameter  beam  pipe  can  supply  a 
magnetic  field  gradient  that  far  exceeds  what  is  needed.  Permanent  magnet  quadrupoles 
can  still  have  adjustable  strengths  depending  on  their  design  [47].  The  total  cost  for  the 
three  quadrupoles  with  the  properties  listed  in  Table  6  would  be  -$40,000  [48].  This 
estimate  includes  the  engineer  design,  the  assembly,  testing,  and  shipping. 

3.  Dipoles 

The  dipoles  were  designed  to  supply  the  transverse  offset  to  the  electron  beam  to 
allow  it  to  be  bent  around  the  mirror  housing  and  also  to  supply  focusing  to  counteract 
the  dispersion  due  to  the  energy  spread.  They  were  designed  to  have  a  nominal  20  cm  arc 
length  leading  to  a  magnetic  field  strength  of -0.056  Tesla.  For  the  first  dipole,  a  rotation 
of  each  end  face  by  15°  is  needed  to  supply  focusing.  For  the  second  dipole  a  rotation  of 
each  end  face  by  -10°  is  needed  to  focus  the  electron  beam  into  the  undulator.  The  total 
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cost  estimate  for  the  two  dipoles  is  -$40,000  with  the  engineer  design,  the  assembly, 
testing,  and  shipping  [48]. 

C.  OPTICAL  CAVITY  BEAMLINE 

The  optical  cavity  beamline  and  adjacent  vacuum  components  must  accommodate 
not  only  the  electron  beam  but  also  the  optical  mode.  The  design  of  this  portion  of  the 
beamline  is  highlighted  in  Figure  50. 


Top  View 


\ 
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Figure  50.  Optical  cavity  beamline  highlighted  in  blue. 


The  optical  cavity  beamline  has  several  main  components.  On  the  two  ends  of  the  cavity 
there  are  two  large  mirror-housing  vessels  that  hold  mirror  mount  assemblies.  These 
assemblies  control  the  longitudinal  /  transverse  positions  and  the  tip  /  tilt  of  the  mirrors  to 
allow  for  proper  alignment.  Directly  on  either  side  of  the  undulator  and  adjacent  to  the 
downstream  mirror  housing  assembly  are  three  crosses  that  allow  for  mechanical 
feedthroughs  to  aid  in  diagnostics  and  alignment.  The  length  of  the  chamber  must  be  such 
that  the  mirrors  can  have  proper  spaeing  to  match  the  optical  pulse  bounce  time  to  the 
eleetron  pulse  repetition  frequency  of  43.75  MHz  as  well  as  being  wide  enough  to  not 
cause  clipping  of  the  optical  mode  in  the  transverse  direction. 

I.  Transverse  Cavity  Dimensions 

The  transverse  size  of  the  components  will  need  to  be  studied  to  guarantee  that  the 
round  trip  losses  stay  below  the  weak  field  gain. 
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a.  Mirror  Size 

The  beam  pipes  surrounding  the  optical  cavity  and  the  mirrors  that  make  up  the 
resonator  must  be  large  enough  to  produce  negligible  clipping  of  the  optical  mode.  The 
large  diffraction  angle  that  occurs  due  to  the  long  wavelengths  of  this  PEL  design 
determines  the  mirror  size  suitable  for  the  Niowave  THz  EEL.  Table  8  shows  the  results 
of  an  analysis  of  the  amount  of  clipping  that  different  wavelengths  will  undergo. 


Mirror 

k  (pm) 

Diameter 

126 

114 

104 

90 

80 

73 

69 

66 

12.7cm 

-0.8% 

-0% 

-0% 

-0% 

-0% 

-0% 

-0% 

-0% 

10.2  cm 

-1% 

-0.2% 

-0% 

-0% 

-0% 

-0% 

7.62  cm 

-1.3% 

-0.5% 

-0.2 

-0% 

Table  8.  Calculated  round  trip  power  loss  due  to  clipping  for  each 
considered  mirror  size  and  wavelength. 


The  clipping  was  calculated  for  three  different  mirror  sizes  that  were  considered 
along  with  a  range  of  wavelengths  that  fall  into  the  Niowave  THz  EEL’s  tunable  regime. 
For  this  analysis,  an  acceptable  round  trip  loss  due  to  clipping  on  the  mirrors  was 
assumed  to  be  ~2%,  as  this  would  allow  a  margin  for  other  areas  of  loss.  The  table  entry 
is  denoted  green  if  the  round  trip  loss  with  the  specific  wavelength  and  optic  diameter  is 
below  that  ~2%  threshold  and  red  if  it  is  not.  A  12.7  cm  (five  inch)  diameter  mirror 
allows  tuning  over  the  entire  considered  spectrum  (~66  pm ->126  pm)  with  less  than  two 
percent  round  trip  losses  from  clipping  on  the  mirrors.  A  10.2  cm  (4”)  diameter  mirror 
size  induces  significant  clipping  of  the  longer  wavelengths,  therefore  limiting  the 
maximum  wavelength  to  about  104  pm.  A  7.62  cm  (three  inch)  diameter  mirror  size 
reduces  this  the  maximum  wavelength  even  further  to  about  80pm.  While  the  12.7  cm 
diameter  option  allows  for  a  full  tunable  range  of  the  EEE,  the  necessary  optomechanical 
components  are  not  as  common  and  go  up  in  price  considerably.  Eor  this  reason,  as  well 
as  for  space  saving  measures,  a  10.2  cm  diameter  mirror  size  is  recommended  for  the 
project. 
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b.  Beam  Pipe  Sizes 

The  beam  pipe  diameter  inside  the  undulator  depends  on  the  neeessary  gap  to 
achieve  a  given  wavelength;  this  analysis  was  discussed  in  Chapter  III.  A  beam  pipe  with 
an  outer  diameter  of  '-'16  mm  and  wall  thickness  of  1  mm  allows  for  wavelengths  --1 00 
pm  and  shorter.  The  beam  pipe  that  encompasses  the  optical  cavity  between  the  undulator 
and  the  mirror  housing  must  be  wide  enough  to  transmit  the  diffracting  optical  mode  but 
should  not  be  larger  than  the  mirrors  as  that  would  not  give  any  advantage.  For  this 
reason  the  beam  pipe  will  be  10.2  cm  in  diameter,  the  same  as  the  chosen  mirror  size. 

2.  Mirror  Mount  Assemblies 

The  mirrors  will  need  to  be  properly  aligned  and  adjusted  to  optimize  the  FEL 
performance. 


a.  Mirror  Mount  Housing  Chamber 

To  hold  the  mirror  mount  setup  on  either  side  of  the  cavity,  a  large  vacuum  vessel 
is  needed  that  supplies  enough  space  for  the  mirror  mount  assembly,  to  allow  a  method 
for  the  light  to  escape,  and  to  allow  for  electrical  feedthroughs  to  control  the  actuators.  To 
meet  these  criteria,  a  large  4-way  cross  will  be  utilized,  as  shown  in  Figure  5 1 . 
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Figure  5 1 .  Mirror  mount  housing  chamber,  from  [49]. 


The  size  chosen  has  a  wall  outer  diameter  of  8”  and  a  centerline  length  (A  in 
Figure  51)  of  7.5”.  They  meet  the  requirements  in  the  following  ways:  they  are  relatively 
inexpensive,  they  have  a  flat  mounting  surface  on  the  bottom  for  the  mirror  mount 
assembly,  they  can  easily  be  attached  to  the  cavity  beam  pipes  with  matching  flange 
styles,  they  have  an  opening  in-line  with  the  cavity  and  the  output  coupler  for  a  window 
to  transmit  the  THz  radiation,  and  they  have  an  opening  to  allow  the  electrical 
feedthroughs  for  the  actuator  controls.  Conflat  style  flanges  that  use  two  knife-edges  and 
copper  gaskets  to  provide  great  seals  that  can  achieve  ultra-high  vacuum  (UFIV), 
typically  greater  than  10'^  Torr.  The  pricing  information  is  listed  in  Table  9. 


Product 

Mirror  housing  chamber 

Company 

Kurt  J.  Lesker  Company 

Model 

C-1000 

Cost/unit 

$2050 

Quantity 

2 

Total 

$4100 

Table  9.  Mirror  housing  chamber  pricing  information. 
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b.  Mirrors 

Due  to  the  lack  of  easily  available  high  reflective  dielectric  coatings  for  the  THz 
regime,  a  metal  mirror  is  recommended.  The  THz  frequency  reflects  well  off  of  metals 
with  a  reflectance  up  to  -99%  (in  power)  for  gold  [50].  Because  of  this,  a  gold-coated 
mirror  is  recommended  with  the  substrate  made  out  of  a  cheaper  material  such  as 
aluminum.  The  final  design  will  be  a  -5  cm  radius  (four  inch  diameter)  gold  coated 
concave  mirror  with  a  radius  of  curvature  of  -1.72  m.  The  middle  of  one  of  the  mirrors 
will  have  a  1  mm  diameter  hole  for  output  coupling.  The  mirror  pricing  information  is 
listed  in  Table  10. 


Product 

Gold  Plated  Mirror 

Company 

Laser  Beam  Products 

Model 

Custom  4”  diameter  mirror,  0.5mm  radius  hole 

Cost/unit 

$680 

Quantity 

2  (one  without  hole) 

Total 

$1360 

Table  10.  Gold  mirror  pricing  information. 


c.  Mirror  Mount  Assembly 

To  support  the  mirror  and  supply  the  necessary  angular  and  translational 
adjustments  the  4”  optic  must  be  held  in  an  optic  mount  assembly.  This  will  consist  of  a 
mirror  mount  to  control  the  angular  alignment  and  a  translation  stage  to  supply 
transverse,  longitudinal,  and  vertical  position  adjustments.  The  components  chosen  must 
be  vacuum  compatible  to  avoid  excessive  outgassing  and  ensure  that  the  vacuum 
requirements  are  met  for  the  whole  system  (-10'^  Torr).  The  stage  must  be  large  enough 
to  supply  the  necessary  length  of  movement  while  not  being  too  large  to  unnecessarily 
increase  the  size  of  the  vacuum  chamber.  To  control  the  mirror  mount  and  translation 
stage  under  vacuum,  motorized  actuators  are  used.  Twelve  actuators  will  be  needed  (three 
per  mirror  mount  and  three  per  translation  stage).  Controllers  will  be  used  to  translate  the 
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orders  from  the  computer  software  to  the  proper  voltages  to  move  the  picomotor 
actuators.  The  parts  price  information  is  listed  in  Tables  1 1-14. 


Product 

Vacuum  rated  4”  mirror  mount 

Company 

Newport  Corporation 

Model 

Custom  version:  U400-AC3K 

Cost/unit 

$2500-$3000 

Quantity 

2 

Total 

$5000-$6000 

Table  1 1 .  Mirror  mount  pricing  information. 


Product 

Vacuum  rated  translation  stage  (12.7  mm  travel) 

Company 

Newport  Corporation 

Model 

9066-XYZ-L-V 

Cost/unit 

$1976 

Quantity 

2 

Total 

$3952 

Table  12.  Translation  stage  pricing  information. 


Product 

Vacuum  rated  picomotor  actuator  (12.7  mm  travel) 

Company 

Newport  Corporation 

Model 

8301-UHV-Kap 

Cost/unit 

$1185 

Quantity 

12 

Total 

$14220 

Table  13.  Picomotor  actuator  pricing  information. 
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Product 

Picomotor  actuator  controller 

Company 

Newport  Corporation 

Model 

8742- 12-kit 

Cost/unit 

$2940 

Quantity 

1 

Total 

$2940 

Table  14.  Picomotor  actuator  controller  pricing  information. 


d.  Output  Window 

The  output  window  material  must  transmit  the  THz  radiation  with  minimum 
absorption  and  be  mounted  on  a  flange  to  be  vacuum  compatible.  A  crystalline  quartz 
window  will  be  utilized  as  it  has  minimum  absorption  in  the  THz  spectrum  and  is  also 
transparent  to  visible  radiation,  which  will  aid  in  alignment.  A  transmission  spectrum  of 
crystalline  quartz  can  be  seen  in  Figure  52. 


Figure  52.  Measured  transmittance  spectrum  for  crystal  quartz  window,  from 

[51]. 
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This  has  a  transmittance  range  of  '-'35%  to  50%  in  the  wavelength  range  that  the  Niowave 
THz  PEL  will  be  operating.  The  noise  in  at  longer  wavelengths  measurement  is  due  to  the 
limitations  of  the  spectrometer  being  utilized.  A  one  ineh  diameter  window  mounted  on  a 
conflat  flange  will  be  used.  The  pricing  information  of  the  window  is  listed  in  Table  15. 


Produet 

Crystalline  quartz  output  window 

Company 

MPF 

Model 

A4977-1-CF 

Cost/unit 

$2350 

Quantity 

1 

Total 

$2350 

Table  15.  Crystalline  quartz  window  prieing  information. 


A  seale  depiction  of  the  final  design  for  the  mirror  housing  can  be  seen  in  Figure 
53. 


Each  element  in  the  mock-up  is  of  the  correct  dimensions  as  the  component  that  will  be 
used.  This  shows  that  the  8”  diameter  mirror-housing  ehamber  is  large  enough  for  the 
transverse  footprint  of  the  translation  stage.  In  Figure  53,  the  optieal  field  will  enter  the 
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chamber  from  the  right  before  refleeting  off  of  the  mirror.  Not  visible  is  the  output 
eoupling  hole,  which  will  allow  the  radiation  to  transfer  through  and  exit  the  output 
window  behind. 

e.  Miscellaneous 

To  transport  the  eleetrieal  signal  from  the  actuator  controllers  to  the  aetuators,  a 
vaeuum  eleetrieal  feedthrough  is  needed.  A  15-pin  feedthrough  will  be  used  because  each 
mirror  mount  assembly  will  need  12  pins  (2  per  actuator)  to  control  the  neeessary 
voltages.  Sueh  a  flange  and  its  neeessary  conneetor  are  seen  in  Figure  54. 


Figure  54.  Example  of  eonflat  flange  mounted  eleetrieal  feed  through  (left)  and 
necessary  eonneetion  (right),  from  [52]  and  [53]. 


These  will  allow  the  movement  of  the  actuators  to  control  the  mirror  alignment.  The 
pricing  information  for  these  parts  is  listed  in  Tables  16-17. 


Product 

Eleetrieal  feedthrough 

Company 

MDC 

Model 

9132017 

Cost/unit 

$376 

Quantity 

2 

Total 

$752 

Table  16.  Eleetrieal  feedthrough  pricing  information. 
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Product 

Electrical  feedthrough  connection 

Company 

MDC 

Model 

9924069 

Cost/unit 

$172 

Quantity 

2 

Total 

$344 

Table  17.  Electrical  feed-through  connection  pricing  information. 

To  aid  in  the  alignment  of  the  optical  cavity  three  diagnostic  ports  will  be  placed 
in  the  optical  cavity  beamline.  They  will  allow  for  three  actuators  to  be  able  to  intercept 
the  cavity.  The  actuators  will  need  precise  control  as  they  must  be  able  to  align  targets  on 
the  ~100  pm  scale  with  the  cavity  axis.  The  option  suggested  for  the  linear  actuator  is 
listed  in  Table  18. 


Product 

Mechanical  linear  feed  through 

Company 

MDC 

Model 

60004-01 

Cost/unit 

$1252 

Quantity 

3 

Total 

$3756 

Table  18.  Mechanical  linear  feedthrough  pricing  information. 


D.  EXIT  BEAMLINE 

After  the  FEE  interaction  in  the  undulator,  the  electron  beam  with  its  induced 
energy  spread  needs  to  be  successfully  guided  into  the  beam  dump.  One  exit  beamline 
geometry  is  highlighted  in  green  in  Figure  55. 
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As  was  stated  earlier,  the  PEL  interaetion  causes  a  large  energy  spread  in  the 
electron  beam  due  to  the  relatively  high  extraction  of  energy.  This  ~20%  energy  spread 
causes  a  large  dispersion  after  the  dipole.  Using  a  combination  magnet  as  before  but  with 
more  extreme  edge  focusing  (larger  canting  of  the  end  faces)  will  be  used  in  an  attempt  to 
counteract  the  large  dispersion  in  angles,  but  it  is  not  very  effective  given  the  magnitude 
of  the  energy  spread.  Therefore,  the  suggested  solution  is  to  place  the  beam  dump  as 
close  as  possible  to  the  last  dipole  (~20  cm),  which  will  keep  the  beam  pipe  and  beam 
dump  a  reasonable  size.  For  this  final  45°,  bend  a  20  cm  long  dipole  with  no  edge 
focusing  a  magnetic  field  of  ~0.1  Tesla  is  needed.  This  will  cost  on  the  order  of  the  other 
dipoles  that  were  quoted  at  -$20,000  [48]. 

E.  ALIGNMENT  TECHNIQUES 

The  alignment  of  the  cavity  mirrors,  electron  beam,  and  undulator  axis  must  be 
sufficient  enough  to  get  some  signal  Ifom  the  EEL.  After  an  initial  signal  is  seen, 
parameters  can  be  further  adjusted  to  maximize  the  signal. 

I.  Transverse  and  Angular  Mirror  Alignment 

The  alignment  system  consists  of  two  targets  on  actuators  that  can  be  placed  into 
the  beam  path,  two  alignment  lasers  with  beam  expanders  and  focusing  optics,  a  pop-in 
mirror  on  an  actuator  to  be  placed  in  the  beam  path,  a  pop-in  mirror  outside  the  cavity  in 
line  with  the  output  coupler,  and  several  cameras  to  view  the  alignment  procedure.  The 
targets  consist  of  triangular  wedges  with  a  1  mm  diameter  hole  drilled  through  each  of 
them.  A  schematic  of  the  system  is  seen  in  Figure  56. 
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Figure  56.  Schematic  of  optical  cavity  alignment  system  for  Niowave  THz  FEL 


First,  prior  to  installing  the  Ml  mirror,  the  undulator  axis  and  the  targets  are 
aligned  using  by  the  undulator  manufacturer,  who  have  mapped  its  magnetic  field.  Once 
the  undulator  axis  and  the  alignment  wedges  are  properly  aligned,  they  can  be  used  as 
targets  for  alignment  of  the  cavity  mirrors.  The  holes  allow  for  an  alignment  laser  to 
propagate  through  the  targets  while  the  angled  sides  allow  the  use  of  digital  cameras  to 
view  any  misalignments.  The  beam  expander  and  focusing  optics  allow  for  the  focus  of 
the  alignment  laser  to  be  shifted,  aiding  in  alignment.  The  steps  are  as  follows. 

1 .  Ensure  that  both  the  upstream  and  downstream  alignment  lasers  (El  and 
E2)  are  collinear  down  the  magnetic  axis  (aligned  with  the  targets  T1  and 
T2)  by  adjusting  the  turning  mirrors. 

2.  Install  the  output  coupler  mirror  (Ml  in  Eigure  56)  and  align  the  hole  with 
laser  beam. 

3.  Removing  P2,  visually  align  the  center  of  M2  along  El . 
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4.  Put  in  P2,  check  alignment  of  L2  with  eavity  axis  and  then  remove  T 1 . 

5.  Using  L2,  view  the  reflection  from  Ml  onto  T2  and  align  the  center  of  the 
reflection  on  the  hole. 

6.  Remove  P2,  insert  Tl,  and  check  alignment  of  LI  with  cavity  axis. 

Remove  T2  and  P2. 

7.  Using  LI  view  the  reflection  from  M2  onto  Tl  and  align  the  center  of  the 
reflection  on  the  hole. 

The  transverse  shift  alignment  of  M2  (step  3)  is  less  stringent  than  other 
toleranees  as  the  mirrors’  spherical  curvature  allow  for  any  small  transverse 
misalignments  to  be  corrected  with  a  tilt,  whieh  is  done  in  step  7.  With  this  teehnique, 
one  ean  achieve  alignment  to  within  -O.!  mm  of  the  cavity  axis  on  the  target  [54].  This 
leads  to  a  cavity  alignment  tolerance  of  -O.!  mrad,  whieh  is  just  inside  the  neeessary 
-0.15  mrad  tolerance  and  will  allow  fine-tuning  adjustments  by  maximizing  the  signal  of 
the  PEL  while  operating.  The  designed  optieal  cavity  beamline  has  three  different 
diagnostic  locations  where  the  two  targets  and  the  pop  in  mirror  can  be  placed. 

2.  Longitudinal  Alignment 

The  EEL  eavity  must  have  the  correct  length  and  desynchronism  on  the  order  of 
the  slippage  length  (-1  mm)  in  order  for  lasing  to  occur.  This  is  fairly  lenient  compared 
to  other  shorter  wavelength  EEE  oseillators.  In  practice  the  oscillator  cavity  length  is 
adjusted  elose  to  the  desired  length  using  simple  measuring  techniques  before  the  cavity 
length  is  seanned,  searehing  for  lasing  operation.  This  is  the  method  that  will  be 
attempted  for  this  EEE  as  adding  other  laser  range  finding  or  survey  methods  would  add 
unneeessary  complexity.  [55]. 

F.  FINAL  DESIGN  IN  3D 

This  design  of  the  optieal  cavity  beamline  and  mirror  housing  units  (see  Figure 
57)  meet  the  requirements  that  were  neeessary  for  the  FEE  to  operate. 
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Figure  57.  Final  design  of  optieal  eavity  beam  vaeuum  struetures 


It  inelude  a  low-loss  optieal  eavity,  the  required  eomponents  for  initial  eavity  alignment, 
and  the  beam  line  eomponents  to  eouple  the  eleetron  beam  into  the  cavity  as  well  as 
ejecting  it  into  the  beam  dump. 
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VI.  CONCLUSION 


Free  eleetron  lasers  are  sourees  of  eoherent  radiation  that  ean  be  designed  and 
tuned  aeross  the  electromagnetie  spectrum.  This  ability  allows  them  to  operate  in  portions 
of  the  spectrum  where  there  are  no  coherent,  powerful  kW-level  sources,  for  example  in 
the  THz  regime.  Due  to  the  nature  of  the  FEL  lasing  medium  (relativistic  electrons)  it  can 
be  scaled  to  high  power.  These  attributes  allow  THz  FELs  to  meet  a  wide  range  of 
applications  including  spectroscopy,  pump-probe  experiments,  and  screening  for  security 
purposes.  THz  EEEs  have  been  created  in  the  past  but  are  generally  large  sized, 
expensive  facilities.  With  the  use  of  a  compact  linac,  a  short  undulator,  and  an  optical 
cavity,  a  relatively  inexpensive,  kW-level,  tunable,  THz  source  can  be  created.  This 
dissertation  presents  and  studies  the  design  of  such  a  THz  EEE. 

The  EEE  was  simulated  using  three  different  codes  to  show  the  effects  of  various 
design  parameters.  These  simulations  show  that  an  electron  beam  with  an  emittance  of 
28mm-mrad  and  energy  spread  of  Jy/y~1.7%  can  be  used  to  create  an  EEE  with  moderate 
weak  field  gain  and  high  extraction.  The  optical  cavity  was  designed  to  operate  over  the 
portion  of  the  spectrum  from  -60  pm  ^126  pm  wavelength. 

The  electron  beam  with  realistic  energy  spread  and  emittance  must  be  transported 
from  the  linac  into  the  undulator  and  must  properly  align  within  the  optical  mode.  To 
control  and  focus  the  electrons,  a  beamline  was  designed  utilizing  quadrupoles  and 
combination  dipole  magnets.  The  combination  dipole  magnets  with  large  amounts  of 
edge  focusing  were  designed  to  overcome  the  dispersion  caused  by  the  electron  beam 
energy  spread  and  to  focus  the  electron  beam  into  the  undulator.  Ultimately,  the  beamline 
components  were  found  to  meet  the  initial  nominal  requirements  of  -0.45  mm  radius  in 
the  center  of  the  undulator  set  out  during  the  EEE  design  and  simulations. 

The  hardware  for  the  Niowave  THz  EEE  was  designed  to  control  the  oscillator 
cavity  mirrors  and  the  electron  beam.  The  optical  cavity  beamline  components  allow  for 
the  mirror  mount  assemblies  to  control  the  mirrors’  angular,  transverse,  and  longitudinal 
adjustments.  Alignment  techniques  were  described,  utilizing  three  different  diagnostic 
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feedthroughs  to  allow  alignment  targets  as  well  as  pop-in  mirrors.  Finally,  the  exit  beam 
line  was  designed  with  a  short  transit  into  a  beam  dump  to  handle  the  ~20%  induced 
energy  spread  generated  by  the  high  extraction  FEL. 

While  this  dissertation  covers  a  large  portion  of  the  major  design  criteria  for  the 
FEL,  several  additional  things  can  also  be  considered  going  forward.  Eringe  field  effects 
must  be  incorporated  into  the  design  of  the  magnetic  field  strengths  of  the  dipoles  and 
quadrupoles.  To  counteract  pulse  elongation  that  occurs  during  the  dipole  offset,  the 
initial  pulse  length  and  the  energy  spread  may  need  to  be  decreased  by  a  factor  of  three  or 
other  bunching  mechanisms  included  into  the  beamline  to  obtain  the  '-Sps  bunch  length  at 
the  undulator  required  by  the  EEL.  Engineering  concerns,  such  as  the  necessary  shielding 
around  the  beam  dump  to  ensure  protection  of  the  surrounding  equipment,  will  need  to  be 
assessed  and  designed.  In  terms  of  the  applications  for  this  EEL,  the  expected  spectrum 
should  be  adjusted  to  meet  the  requirements  set  out  by  the  application. 

This  design  produces  an  ~2  kW-level,  compact,  THz  EEL  for  a  wide  range  of 
applications.  The  compact  size  and  relatively  inexpensive  components  open  up  the  use  of 
THz  radiation  to  many  new  users. 
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